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ABSTRACT 


This report presents the results of several 
analytical studies related to Space Shuttle Navigation. 
It is divided into two Volumes, the first dealing 
with studies related to the addition of NAVSTAR Glo- 
bal Positioning System user equipment to the Shuttle 
avionics suite, and the second dealing with studies 
of the baseline avionics suite without GPS. The GPS 
studies center about navigation accuracy covariance 
analyses for both developmental and operational phases 
of GPS as well as for various Orbiter mission phases. 
The baseline navigation system studies include a 
covariance analysis of the Inertial Measurement Unit 
calibration and alignment procedures, postflight 1MU 
error recovery for the Approach and Landing Phases, 
on-orbit calibration of I MU instrument biases, and a 
covariance analysis of entry and prelaunch naviga- 
tion system performance. 
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1 . 


INTRODUCTION 


This is Volume II of a two volume final report on 
Space Shuttle Navigation Analyses. Volume 1 is concerned with 
the addition of GPS user equipment to the Shuttle avionics 
suite, while this volume treats several topics related to the 
baseline navigation system for the Orbiter. An introduction 
to the entire report, relating the current studies to earlier 
work, is to be found in Volume I. The remainder of this intro- 
duction, extracted from that in Volume I, refers only to the 
work reported here. 


1.1 SHUTTLE BASELINE NAVIGATION SUITE 

Figure 1.1-1 illustrates the various mission phases 
of a Space Shuttle flight and helps tie together the various 
studies of the baseline navigation system reported on in this 
volume. The elements of the baseline navigation system include 
Inertial Measurement Units (IMU's) with their inertial sensors, 
barometric altimeters, TACANS, and a Microwave Landing System 
(MLS). The MLS is not addressed in the present work; it has 
been studied in prior work. The other elements of the base- 
line navigation system are all involved in the studies of this 
volume . 


1.2 SHUTTLE BASELINE NAVIGATION SUITE 

There are four areas in which performance of the baseline shuttle 
navigation system was addressed in this study. These areas 
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ORBITAL COAST 

ACCUCflOUWIIt AND 
GYROSCOPE BIAS CALIBRATION 





Figure 1.1-1 Baseline Navigation Studies Overview 

are related in that they all focus on the performance charac- 
teristics of the Orbiter Inertial Measurement Units (IMU's). 
Three of them involve calibration of the IMU, while the fourth 
centers on navigation performance. 

Preflight calibration and alignment of the IMU's has 
been studied in prior contracts (Ref. 1), but this report pro- 
vides the first covariance analysis of the current Orbital 
Flight Test (OFT) cal/align procedure (covariance analysis 
results for two prior versions of the OFT cal/align procedure 
were also generated under this contract -- see Refs. U and 5). 
The previous efforts utilized monte carlo techniques to gen- 
erate an overall performance projection for the Approach and 
Landing Tests (ALT), and to identify potential risk areas. 

The current effort, reported in Chapter 2, uses covariance 
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techniques to generate overall performance projections, to 
develop a detailed error budget for the OFT missions, and to 
identify the major error mechanisms that limit calibration and 
alignment performance. 

ALT provided NASA's first opportunity to evaluate 
both the navigation performance achievable with the 1MU and 
the true calibration and alignment uncertainties in a flight 
environment. The study of postflight IMU error recovery re- 
ported in Chapter 3 provides a preliminary assessment of the 
potential accuracy of recovering IMU errors from ALT. The 
performance specification for the Shuttle IMUs requires that 
they maintain their prelaunch calibration acucracy for a peri- 
od of approximately 15 hrs. After lengthy periods on-orbit, 
the IMU gyros and accelerometers will need recalibration of 
their observable characteristics (primarily biases) in the 
weightless environment. In Chapter 4, the best method for 
performing such a calibration is addressed, and procedures, 
timings, and analytical techniques are recommended based on 
laboratory models of the instruments. 

In Chapter 5, an assessment of navigation performance 
during entry but prior to the landing phase is provided. This 
work is a follow-on to work completed under the previous con- 
tract (Ref. 4), but is based on more recent, optimistic assum- 
ptions of initial errors and sensor performance, and a revised 
measurement schedule for the baro altimeter. A candidate, 
baseline navigation filter is examined in detail. 

Chapter 6 is a summary chapter, highlighting the con- 
clusions from all of these baseline navigation system studies. 
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2. IMU CALIBRATION AND ALIGNMENT COVARIANCE ANALYSIS 


The performance analysis for the Space Shuttle Inertial 
Measurement Unit (IMU) calibration and alignment software has 
been conducted in two phases. The first phase, completed under 
a previous contract period (Kef. 1), used monte carlo techniques 
to generate an overall ' erformance projection for the Approach 
and Landing Test (ALT) calibration and alignment algorithm 
(Ref. 2) and to identify potential risk areas. The second 
phase uses covariance analysis techniques to develop a detailed 
error budget for the Orbital Flight Test (OFT) algorithm (Ref. 

3) and to identify the major error mechanisms which limit cali- 
bration and alignment performance. For those error parameters 
for which the OFT algorithm cannot provide the required calibra- 
tion accuracy, the results of this phase provide a clear indi- 
cation as to the algorithm modifications which should be made 
in order to achieve the desired performance. 

An error budget showing the contribution of all error 
sources to calibration and alignment performance for an earlier 
version of the cal/align algorithm was presented in Refs k and 
5. That error budget was applicable to the OFT algorithm de- 
fined by Ref. 3. dated December 1976. The cal/align error 
analysis discussed in this report incorporates all changes 
made to the OFT algorithm through December of 1978 wnich may 
affect cal/align performance. The most important of these 
changes are: 


• New 13-position Hangar Cal A and 7- 

position Preflight Cal A sequences de- 
signed to avoid the effects of the tran- 
sients induced when switching accelero- 
meters between high and low’ gain (Ref. 6) 
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• New Hangar Cal A and Pre flight Cal A 
calibration equations to accommodate the 
13-position and 7-position sequences 
(Ref. 6) 

• Modified Hangar Cal A calibration equa- 
tions for the high gain accelerometer 
asymmetry errors (Ref. 7) 

• Changes to the cal/align sequences for 
Hangar Coarse Alignment, Preflight Coarse 
Alignment, and Gyrocompassing 

• Incorporation of the accelerometers' 
double dead band characteristics and its 
software "fix.” An evaluation of the 
impact of LSF filter estimation accuracy 
of this effect was reported in Ref. 8 

• An update of the 1MU error model data 
base based upon inputs from Singer- 
Kearfott Division (SKD) and Rockwell 
International (RI) (Ref. 13-19, and 30) 


An overview of the calibration and alignment mechani- 
zation for OFT is presented in Section 2.1. Section 2.2 con- 
tains a general description of the mathematical structure used 
in the covariance simulations, and Section 2.3 summarizes the 
truth model states and error sources used in the covariance 
analysis. An error budget showing the contribution of all 
error sources to calibration and alignment performance, as well 
as a discussion of the results obtained, is presented in Section 
2.4. 


2.1 CALIBRATION AND ALIGNMENT OVERVIEW 

The calibration and alignment of the Space Shuttle 
IMUs consists of two parts: Hangar Calibration, and Preflight 

Calibration and Alignment. Hanger Calibration will be performed 
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4 ^ ^ 

in a low-vibration environment 1 within 5 to 14 days prior to 
launch, and will have a total duration of 9.6 hr. Preflight 
Calibration and Alignment will be initiated at the launch pad 

-jfr " n_l 

within 3 to 15 hr before launch, and will have a total dura- 
tion of 3.1 hr. The calibration and alignment schedule during 
both the hangar and the preflight mission phases is illustrated 
in Fig. 2.1-1. 

The cal/align mechanization consists of steering the 
platform to a series of discrete positions at which torquing 
rates are applied to compensate for earth rate, and in some 
cases "excess" rates are applied to torque the platform through 
a prescribed trajectory relative to a local vertical (North, 

West, Up) coordinate system. At each of these positions, the 
accelerometer outputs and/or the resolver outputs are sampled 
and a least squares fit (LSF) filter is used to estimate the 
acceleration along the up axis, the platform tilts and tilt 
rates about the north and west axes, and/or the platform atti- 
tude and attitude rate. Once all data (vertical accelerations, 
tilts, etc. for each discrete position) have been collected, a 
number of IMU errors are estimated and the corresponding error 
compensation parameters are updated. Between discrete positions, 
the platform is slewed at a high rate and no data is collected 
during this motion. 

The cal/align mechanization is comprised of six diff- 
erent sequences. Each of these consists of several positions 
and is used for the calibration, alignment, or verifciation of 
several IMU errors as outlined in Table 2.1-1. 


For the OFT mission, Hangar Calibration is performed on the 
launch pad, this is reflected in the results reported here. 

•!« 

These times are approximate and subject to modification. 
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Figure 2.1-1 Calibration and Alignment Schedule 
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2.2 COVARIANCE ANALYSIS PROGRAM DESCRIPTION 
2.2.1 Program Overview 

The structure of the covariance analysis program is 
the same as that of the OFT IMU software. The program is 
organized into seven ' eparate sequences with an executive pro- 
gram which defines the interface between those sequences (Fig. 
2.2-1). The input to the Coarse Alignment sequence is a data 
file which defines the IMU status (values for all error co- 
efficients, gyro torquing rates*, etc.), a priori calibration 
and alignment estimates, and the transformation matrices be- 
tween the platform coordinate system and the vehicle-, earth- 
and inertially-fixed coordinate systems. The output of the 
Coarse Alignment sequence is a new data file updated by the 
Coarse Alignment estimates . This file is then the input to 
Hangar Cal A and so on. 

The program preserves the correlation between esti- 
mation errors in different sequences and has the capability of 
simulating a complete calibration and alignment procedure 





Figure 2.2-1 Structure of Covariance Simulation 

for Performance Analysis 
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from beginning to end. Alternatively, the initial condition 
data for any sequence can be used (and modified if desired) 
for performing special studies for any other given sequence. 

2.2.2 State and Covariance Equations 

The LSF filter used in the cal/align software is a 
recursive filter, i.e., the tilt and tilt rate estimates gener- 
ated at a particular time are a function only of the current 
acceleromter outputs and the estimates generated at the previous 
time. This recursive property permits the estimation procedure 
to be modeled as a (suboptimal) state estimator for the IMU 
error parameters based on a linearized model of the system 
dynamics. Because of the fact that the IMU paramter estimates 
are updated only at the end of each cal/align sequence, how- 
ever, it is necessary to introduce temporary storage states 
into the filter equations. The form of these equations and 
the approach used in the program to generate the accompanying 
covariance equations are discussed in this section. 


The state and covariance error equations are summarized 
in Table 2.2-1. Equation 2.2-1, found in the table, models 
the error propagation for the linearized IMU system. This 
equation represents the fundamental operational mode of the 
IMU including its associated software during the Hangar and 
Preflight phases (the Ground Sequence, Ref. 3). The state 
vector x is composed of three parts: 


Xj - the least squares fit (LSF) filter states 

used to estimate the platform accelerations 
(tilts) and acceleration rates (tilt rates) 
by measuring the platform velocity outputs, 
and/or to estimate the platform attitude 
and attitude rate by measuring the gimbal 
angles . 
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TABLE 2.2-1 

SYSTEM AND COVARIANCE ERROR EQUATION SUMMARY 


CAL ALIGN STEP 


STATE EQUATION 


ESTIMATE TILTS AND DRIFTS 
PROPAGATE !MU ERRORS 


TORQUE OUT TILTS 
RESET FILTER 


STORE ESTIMATES 
RESET FILTER 


UPDATE I MU PARAMETERS 
RESET STORAGE STATES 



* \v 

* 2 . 2 -n 




COVARIANCE EQUATION 


F P - P F T - Q (2.2-2' 




\n^re :r.e vectors ar* deiined bv 

LSF Fi Per States 

\ * \ 0 * Storage Sta:*.« 

I MU Error States 
** m 

a:td the matrices are j:\en bv 


-1 

u » 0 


Measurement Noise 


Process Xoi*e 


i.. 

;i. iUNa 



0 0 

KH 

F « 

A 

0 


0 0 



ERROR COVARIANCE 


? P P 

'll *12 *13 

P * r p p 

12 r 22 23 


PLATFORM LEVELING 


0 0 0 

L ® 0 1 3 

-L, 0 I 


p,~ P-- 

1 j «. j j 3 

DATA STORAGE 


0 0 0 

c„ I c„ 

0 0 l 


NOISE COVARIANCE 


R i: 0 0 

Q * o o o 

0 0 Q„, 


PARAMETER UPDATE 


A 

0 

« 

c 

0 


0 

0 

» 


I 


x 9 - the storage states, used to store some of 
the filter estimates, x, , which will be 
used to update 1MU error paramters at a 
future time. 

x^ - the system states, including platform ve- 
locity errors, misalignments, and all I MU 
error states being considered in the present 
cal /align evaluation. 
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Equation 2.2-1 is driven by white process noise represented by 
the vector w. The components of w are resolver measurement 
error and the IMU process noise w^. There are no errors 
associated with the transfer of filter states x^ into storage 
states * 2 ' The propagation of the statistics associated with 
the state vector x is governed by Eq. 2.2-2. The details about 
the elements in and H are given in Appendix B of Ref. 9. 
Details about the remaining matrices are given in Ref. 10. 


At the end of each estimation period some of the fil- 
ter estimates x^ are stored in the temporary memory states 
and. simultaneously, the filter states are reset to zero. This 
storage and reset is represented by Eqs. 2.2-5 and 2.2-6. 
Depending on the particular point in the cal/align process, 
the next step .nay be one of the following: 

• Update one or more of the 1MU errors, 
x^, using one or more of the estimates 
stored in x«» and reset the storage states 
in x^, whicn are not required for future 
1MU Srror updates. This step is repre- 
sented by Eqs. 2.2-7 and 2.2-8 

• Rotate the platform about the horizontal 
axes by an amount equal and opposite to 
the filter estimates of the platform 
misalignments from level. This step is 
represented by Eqs. 2.2-3 and 2.2-9 

• Initiate the next estimation of platform 
tilts and drifts and IMU error propaga- 
tion, as described above and represented 
by Eqs. 2.2-1 and 2.2-2. 

The cal/align evaluation performed by TASC is based 
on computation of the covariance Eqs. 2.2-2, 2.2-9, 2.2-6 and 
2.2-8. Since the corresponding state Eqs. 2.2-1, 2.2-3, 2.2-5 
and 2.2-7 are all linear, it is possible to initialize the 
matrices P and Q with only those error sources of interest 
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at a particular time and to obtain, at the end of the run, the 
sensitivity of the cal/align mechanization to that particular 
error source or group of errors. 


2.3 TRUTH MODEL DESCRIPTION 

Generation of an error budget for the cal/align 
software requires that a "truth model" describing the real 
world error sources for the cal/align process be defined. The 
description must include a complete list of states and error 
sources, and a data base. The truth model used to evaluate 
the Space Shuttle IMU cal/align performance is presented in 
this section. This model incorporates all changes to the 
baseline OFT algorithm mentioned at the outset of this Chapter. 

2.3.1 States and Error Sources 

The truth model states and error sources used in 
evaluating the IMU cal/align mechanization are listed in Table 
2.3-1, which divides them into three major categories: 


• 

Category 

1 

Platform misalignments and 
calibrated IMU error states 

• 

Category 

II 

Noncal ibrated IMU error 

states 

» 

Category 

III 

Random and quantization 

errors 


These three categories contain the error sources 
which will be used to generate the baseline error budget . The 
first category corresponds to the platform misalignments 


The "truth model" is a mathematical model of all potentially 
significant error sources and the way they affect the cal/ 
align performance in the real world. 
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TABLE 2.3-1 

1MU TRUTH MODEL STATES AND ERROR SOURCES 


T-004 la 


ERROR SOURCE SAME 
AND 

CROUP NUMBER 

NUMBER 

OF 

STATES 

NUMBER 
OF ERROR 
SOURCES 

ERROR 

MNEMONIC 

I. PLATTORM MISALIGNMENTS 

AND CALIBRATED ERRORS 




1. 

Platform Misalignment 

3 

0 


2. 

Gyro Bias Drift 

3 

3 

DFX.Y.Z 

3. 

Gyro Mass Unbalance 

5 

5 

DIX .2 . DSX.2 .DOZ 

4. 

Gyro Torquer Scale Factor 

3 

3 

KTX.Y.2 

5. 

Gyro Misalignment 

6 

6 

BXY.XZ ZY 

6. 

Accel. Bias -Low Gain 

3 

3 

KOX.Y.Z 

7 . 

Accel. Scale Factor -Low Gain 

3 

3 

K1X.Y.Z 

8. 

Accel. Bias -High Gain 

3 

3 

KOHX.Y.Z 

2. 

Accel. Scale Factor -High Gain 

3 

3 

K1HX.Y.Z 

10. 

Accel. Scale Factor Asvmmetry 

6 

6 

KSX.Y.Z.KSHX.Y.Z 

11. 

Accel. Nonorthogonality 

3 

3 

DEL7YX , ZX . ZY 

12. 

Accel, and Gyro Misalignment 

2 

2 

MXZ.MYZ 

13. 

Resolver Offset. Gimbal 
Nonorthozonalitv 


6 

IRO.PO.AZO.DP.a.b 

1 II. NOXC AL I SRATED IMU ERRORS 




14. 

Gyro Anisoelasticities 

SSI 

11 

DI2X.Z.DS2X. Y.Z 
D ISX . Z . DOSX , Z , D IOX . Z 

15. 

Gyro Output Axis Mass Unbalance 

2 

2 

DOX.Y 

16. 

Gyro Heading Sensitive Drift 

IS 

18 

HS ( X . V t Z ) ( OR . P , AZ ) 

1 1 . 

Gyro Thermal Transient Drift 

3 

2 

TTX.Y,Z 

IS. 

Accel. Nonlinearity - 2 nd Order 

6 

6 

K2X.Y.Z.XZ.YZ.ZX 

19. 

Accel. Nonlinearity - 3 rd Order 

3 

3 

K3X.Y.Z 

20. 

Accel. Attitude Sensitive Bias 

18 

18 

KOAS(X.Y.Z) 

(OR.P.PZ) 

21. 

Accel. Attitude Sensitive Scale Factor 

IS 

18 

K1AS(X.YZ ) 
(OR.P.AZ) 

22 . 
23. 

Outer Roll Offset and Misalignment 
Resolver Harmonic 

V 

24 

3 

24 

0R0.-Y.1Z 
1 st . 4 th . 8 th . 9 th 

! II . RANDOM AND QUANTIZATION ERRORS j 




24. 

Gyro Randomness and Quantization 

O 

6 


25. 

Accel. Randomness and Quantization 


6 


26. 

Resolver Randomness and Quantization 


7 


27. 

Vehicle Motion 

7 

2 


TOTALS 

B 

173 



'Only thre* errors are calibrated. 
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and those IMU error sources which are calibrated during one or 
more calibration sequences. Category II corresponds to addi- 
tional IMU error sources which are not estimated or calibrated 
during either hangar or preflight calibration. Category 111 
contains those states which are driven by random errors, 
including instrument randomness, quantization effects, and 
vehicle motions induced by wind gusts. 

All error source categories are divided into smaller 
sets, each of which is associated with a group number. These 
group numbers are useful in defining the detailed truth model 
equations in Appendix B of Ref. 9. In addition, each line of 
the error budget generated in this study will correspond to 
the contribution from all error sources in one particular group 
rather than each error taken one at a time. This simplifies 
the error budget table and places the various error source 
groups in better perspective. 

2.3.2 Truth Model Data Base 

Generation of a detailed error budget requires the 
numerical values for all truth model error sources. Two sets 
of errors must be specified: the initial calibration error 

(instrument error minus software compensation term) at the 
beginning of the hangar calibration phase, and the changes 
of the instrument errors during the time elapsed between the 
hangar calibration and the preflight calibration due to the 
turn-off turn-on errors. This data is summarized in Tables 
2.3-2 through 2.3-4; all error terms which were modified from 
those used by TASC in an earlier cal/align performance evalua- 
tion (Refs. 4 and 5) are marked with a "t . ” 
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TABLE 2.3-2 

1MU TRUTH MODEL DATA BASE FOR CALIBRATED ERRORS 


T-3276 



KliKOK SOtINCK NAMh AND 
GKOUM NUMBER 


SIAN HAND DEVIATION AND SOURCE OF DATA 



UNITS 

A'l BK< 
OK IIAI 

".INNING 
MtiAK GAI. 

THUN 

TUKN 

INSTAK1 

-OFF 

-ON 


> t 

Gym Idas Drift 

sec/sec 

0.035 

Net . 

J 

0.0»5 t 

BBSS 

Kef 

11 

i. 

Gyro Ma^s Unbalance 
Input Axis Acceleration 

SeC/bec/g 

O.'tOO 

Net . 

1 

0.0251 

Ref 

n 


Spin Axis A« celt- rat ion 

sec/stc/g 

0 . 020 

Kef . 

1 

0.025! 

0.025' 

Kel 

11 


Oul|*ul Axis AtielfiMl ion 

see/*ec/g 

0.025 

Kel . 

1 

Kel 

11 

/». 

Gyro Torque r Scale Factor 

ppn. 

/.oo 

Kef. 

1 

100 

Kef 

1 

f >. 

Gyro Mi sa 1 igniuent 

ST? 

120 

Ref. 

1 

oo 1 

Kef 

u 

It. 

Accel. Bias - l.t»w Gain: 
X.Y 

Mg 

100 

Kef. 

1 

SO 

Ket 

n 


L 

MB 

200 

Kel . 

1 

50 

Kef 

. 11 

/. 

Accel. Si*. lie Factor - l.ow Gam 

ppm 

200 

Kef. 

1 

100 

Ket 

. 11 

K. 

Accel, Nias - High !«ai »: 
X.Y 

Mg 

100 

Kef. 

1 

50 

Kel 

. M 


/ 

Mft 

200 

Kel . 

1 

f >0 

Kel 

. 11 

T. 

Accel. Scale factor * High Gain 

ppm 

200 

Kel . 

1 

100 

Kef 

. 11 

10. 

Arcd . Scale factor Asyniwet ry : 
High and l.ow (hi in. Gorrelaled >. 

ppm 

200. 
S V 

Ket. 

1 

20 

See 

Text 


High and Low Gain, Uneorre 1 at ed 

ppm 

Sec Text 

See 

Text 

t 1 . 

At cc 1 . Nonort liogona lily 

t iei* 

00 

Kef. 

1 

21 T 

Ket 

. If 

12- 

Act el. and Gy i o M i sa l i ghincn l 

sec 

00 

Kel . 

1 

15 1 

Kel 

. 11 

i ;. 

Ki'hU 1 V*e r Ol 1 set 

sec 

200 

Kef . 

1 

20 

Ke 1 

. 1 


G i mini 1 N tutor I Logon a I 1 1 y 

sec 

00 

Kel . 

1 

0 

Kel 

. 1 


•ill. high low ^iiii mil liiiWi'Ui' m«iU' 1 act or difference in not calibrated. 


lititse irioi terms dillcr 1 1 oiii those kivcii in He! . 4; they arc discussed in the text. 


The numerical values described above are required to 
construct F, P and Q in Eq. 2.2-2. The specified elements of 
F, P and Q which must be assigned are: 

the initial 1MU error covariance matrix 

the IMU turn-off turn-on instability 
error covariance matrix 

the process noise covariance matrix 


*33 

AP 33 


<33 
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TABLE 

1MU TRUTH MODEL DATA BASE 


2 . 3-3 

FOR NONCALI BRATED ERRORS 


T-327? 


ERROR SOURCE NAME 

VNITS 

STANDARD 

DEVIATION 

DATA 

SOURCE 

1*. 

Gyro Anisoelasticity : 





Input- Input Anisoel. Drift 

sec/sec/g 2 

0.005 

Ref. 12 


Spin-Spin Anisoel. Drift 

sec/sec/g* 

0.005 

Ref. 12 


Input-Spin Anisoel. Drift 

Sec/sec/g 2 

0.025 

Ref. 12 


Output-Spin Anisoel. Drift 

sefc/sec/g 2 

0.025 

Ref. 12 


Input-Output Anisoel. Drift 

secysec/g 2 

0.005 

Ref. 12 

i5 . 

Gyro Output Axis Mass Unbalance 

sec/sec/g 

0.005 

Ref. 12 

lb. 

Gyro Heading Sensitive Drift 

sei/sec 

See Table 2.3-4." 

See Text 

17. 

Gyro Thermal Transient Drift 

sec/sec 

0.01 e* t/12 ° 11 " 

See Text 

IS. 

Accel. Nonlinearity - Second 
Order 


Low Gain High Gain 



Quadratic 

MR/ g^ 

15 I50 t 

See Text 


Product Nonlinearity 

MR/g‘ 

10 100* 

See Text 

19. 

Accel. Nonlinearity - 
Third Order 

3 

M«/« 

5 50* 

See Text 

20. 

Accel. Attitude Sensitive Bias 

Mg 

See Table 2.3-^ 

' 

See Text 

21. 

Accel. Attitude Sensitive 
Scale factor 

1 

ppm 

See Table 2.3-A* 

See Text 

22 . 

Outer Roll Offset and 


200 

t 

Ref. 1 


Misalignment 

sec 

23. 

Resolver Harmonic 


7.6 (2 > 



i fiC Harmonic 

tel 

Ref. 12 


8 th Harmonic 

sec 

19.0 (2) 

Ref. 12 


9 th Harmonic 

sefc 

4.2 <2 ' 

Ref. 12 


1 6 th Harmonic 

sec 

20. 0 t2 ’ 

| 

P.ef . 12 

lu. 

Gyro Randomness 

sec/sec*^ 

2. 69x10* 5 

See Text 

25 . 

Accel. Output Noise 

Mg | 

5.0 

Ref. 1 


Velocity Quantization- Low Gain 

ft/sec 

See Text 

See Text 


Velocity Quantization -High Gain 

ft/sec 

See Tex. 

See Text 

26. 

Resolver Output Noise 

sec 

: 

See Text 


Resolver Quantization 

SSt 

5.8 

See Text 

27 1 

Vehicle Motion 


See Text 

See Text 


il> t « seconds after completing a cluster repositioning. 

(2> Standard deviation of maximum amplitude. 

’These error terms differ from those given ir. Ref. *•: they are discussed m *.h. text. 
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TABLE 2.3-4 

IMU TRUTH MODEL DATA BASE FOR HEADING SENSITIVE 
GYRO AND ATTITUDE SENSITIVE ACCELEROMETER ERRORS 


T-3278 





STANDARD DEVIATION 1 

GYRO OR 


ERROR 

OF ERROR COEFFICIENT I 


rtuD ki 




AUL£.i. . 

u 1'IDAL 

TERM* 




AX I S 


16. Gyro 


21. Accel. 




Drift 
( ’see'/ sec) 

HU 

Scale 

Factor 




(Mg) 

(.ppm) 


Roll 

sin (2<j>) 

0.004 

4.3 



cos (2o) 

0.004 

4.3 


V 

Pitch 

sin (2d ) 

0.004 

4.3 


A 

cos (26) 

0.004 

4.3 

ifflBI 


Azimuth 

sin v2u>) 

0.008 

3.7 



cos (2u) 

0.008 

3.7 



Roll 

sin (2o) 

0.004 

1.1 

3.4 


cos (2o) 

0.004 

1.1 

2.4 

Y 

Pitch 

sin (26) 
cos (26) 

oo 

oo 

oo 

2.2 

2.2 

6.7 

6.7 


Azimuth 

sin (2o) 
cos (2u) 

0.008 

0.008 

3.7 

3.7 

11.5 

11.5 


Roll 

sin (2o) 

0.008 

2.2 

6.7 


cos (2o) 

0.008 

2 ! 2 

6.7 


Pitch 

sin (26) 
cos (26) 

0.008 

0.008 

4.3 

4.3 

1 

' 5:i : 


Azimuth 

sin (2u) 
cos (2v> 

0.006* 

0.006 t 

3.5 

3.5 



, 6. ana .• represent roll, pitch, and azimuth gimbal angle, respectively. 

T These error terms differ from those given in Ref. 4; they are discussed 
in the text. 


" the IMU truth model error dynamics matrix 
Rj^ - the measurement error covariance matrix 

Elements of P33* Q33» anc * P33 are normally chosen together to 
define random processes with desired properties. For instance, 
a first-order markov process (exponentially correlated random 
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process) is modeled with a specified rms value, o^ , and corre- 


lation time, x^. In 
elements are: 

this case the relevant 

truth model matrix 

P 33 . . 
JJ i i 

= a? 
1 

(2.3-1) 

Hi 

* 

(2.3-2) 

933 u 

* 2o H 

(2.3-3) 


A constant error source can be considered a special case of 

~ Qo o = 0 . 


the above with = * and 


33 

ii ii 


Special characteristics of some of the error groups 
shown in Tables 2.3-2 through 2.3-4 are discussed below. 

Ail error terms whose numerical values were modified from 
those used by TASC in an earlier cal/align performance evalu- 
ation (Refs. 4 and 5) are pointed out; these are marked with 
a "+" in Tables 2.3-2 through 2.3-4. 

Preliminary versions of Tables 2.3-2 through 2.3-4 
were reviewed with SKD (Ref. 13). Based on SKD's recommendations, 
the acclerometer scale factor pitch gimbal sensitivity terms 
in Tables 2.3-2 through 2.3-3 were reduced; these are discussed 
under the heading "Groups 16, 20, and 21". 

Group 10; Acceleromter Scale Factor Asymmetries 
(Table 2.3-2) - The model for Group 10 in Ref. 4 assumed 
uncorrelated low and high gain scale factor asymmetry errors. 

Since both types of errors were being calibrated, the magni- 
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tude and degree of correlation of these errors was not impor- 
tant. More recent test data, however, has indicated the 
following (Ref. 14): 

e The high and low gain asymmetry errors 
are partially correlated 

e The difference between high and low gain 
asymmetry errors is at most 100 ppm (3o) 

e Certain accelerometer errors (including 
Croups 18 and 19 •• accelerometer non- 
linearities) are ten times larger in 
high gain than in low gain, which can 
lead to difficulties in calibrating the 
high gain scale factor asymmetry. 

A Software Change Request (Ref. 7) was formulated to remedy 
the latter problem by using only low gain measurements to cali- 
brate both high and low gain accelerometer asymmetry errors. 
Accordingly, the new model for Group 10 in Table 2.3-2 exhibits 
a 200 ppm rms correlated error with a 20 ppm turn-off turn-on 
instability for the low gain asymmetry error. The high gain 
asymmetry error equals the low gain error plus an uncorrelated 
error of 33 ppm (100 ppm, 3o) with no turn-off turn-on insta- 
bility. Note that now the uncorrelated part of the high gain 
asymmetry error is not calibrated . This shortcoming is deemed 
less serious than the errors that are introduced when attempt- 
ing to calibrate the high gain asymmetry errors using high 
gain accelerometer data (Ref. 14). 

Turn-off Turn-on Instabilities (Table 2.3-2) - All 
error terms marked with a "t” were increased in magnitude from 
those used in Ref. 4. The new values are the long-term sta- 
bility errors specified in Ref. 8. These changes were recom- 
mended by SKD (Ref. 15). The values for error Group 10 are 
discussed in the foregoing paragraph. 
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Groups 16. 20, and 21: Heading Sensitive Drift and 
Attitude Sensitive Bias and Scale Factor (Table 2.3*4) - Each 
of these INU error sources is modeled as the sum of three, 
two-cycle sinusoidal errors, one each for the outer roll, pitch 
and azimuth gimbal angles. For the baseline error budget pre- 
sented in Sections 2.4.1 and 2.4.2, the phase of each sinusoid 
is assumed random, and thus each sinusoid is modeled as the 
sum of a sine term and a cosine term that are independent and 
have equal rms magnitudes as listed in Table 2.3-4. 


The majority of the error magnitudes shown in Table 
2.3-4 were deduced from SKD test data reported in Ref. 16. 
Based on several discussions with R1 and SKD, several errors 
were modified as follows: 


• The rms Z-gyro azimuth heading sensitivi- 
ties were reduced from 0.025 to 0.006 
sec/sec. This change results from an 
improved heater design. It is based 
upon test data reported in Ref. 30. 

• The Z-acceleromter scale factor and bias 
azimuth attitude sensitivities were re- 
duced by a factor of 2, to 11.0 ppm and 
3.5 pg, respectively. These new values 
are based on test data reported in Ref. 

30. Nevertheless, these errors continue 
to be significant contributors to cal / 
align inaccuracies because the Imposition 
calibration sequence is Highly sensitive 
to azimuth attitude sensitive errors^ 


• The rms X- and Z-accelerometer scale 

factor pitch attitude sensitivity errors 
were reduced to 7.1 ppm (i.e., total 
magnitude of 10 ppm). This change arises 
because, according to Ref. 7, these error 
terms are now dynamically compensated by 
the IMU software. According to R1 (Ref. 
14) the rms residual errors of these 
terms are less than 10 ppm. 
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Group 17: Gyro Thermal Transient Drift (Table 2.3-3) - 

A new magnitude and time constant were recommended by SKD 
(Ref. 17) for the gyro thermal transient drift error. Because 
of the smaller time constant (originally 600 sec), the effect 
of this error group will be diminished. 

Groups 18 and 19: Acceleromter Nonlinearities 

(Table 2.3-3) - SKD has pointed out (Ref. 17) that the second- 
order and third-order accelerometer nonlinearities are ten 
times larger in high gain than in low gain. With this change 
these error groups will become significant contributors to 
calibration inaccuracies . 

Group 24; Gyro Randomness (Table 2.3-3) - The gyro 
randomness error model used by TASC assumes a random walk (Ref. 
10). The growth rate has been adjusted to yield the same error 
growth over periods of up to ten minutes (the maximum LSF filter 
period) as the SKD markov model presented at the February 1977 
IMU Verification Meeting (Ref. 18) and is equivalent to a vari- 
ance of 3 x 10* (deg/hr r/hr. (The ramp model presented in 
Ref. 18 is insignificant over periods of up to one hour, relative 
to attitude sensitive and thermal transient drifts.) 

More recently, SKD has developed a new gyro randomness 
model based on IMU test data (Ref. 19). Preliminary calcula- 
tions performed by TASC indicate that the average LSF filter 
estimation errors for filtering times of up to 50 minutes (the 
maximum estimation period) with the SKD model would exceed 
those with the current TASC model (Ref. 4) by some 22%. Since 
additional information on the SKD model has not yet been obtained, 
the TASC model has been retained. 
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Group 25: Accelerometer Output Noise (Table 2.3-3) 

The accelerometer output noise model consists of a markov 
process and a quantization error. The markov process has an 
rms error level of 5 p£ and a correlation time of 0.4 sec 
(Ref. 1). The velocity quantization error with the acceler- 
ometers operating in the low gain mode is determined as 
follows: 


One AV pulse * 0.0344 fps (2.3-4) 

Each time an element of the final velocity vector is sampled, 
quantization contributes an error ranging from -0.0172 fps to 
+0.0172 fps. Assuming that this total velocity error is inde- 
pendent from velocity sample to velocity sample, each AV pulse 
count has a random error selected from a triangular distribu- 

i 

tion shown in Fig. 2.3-1. Although the total velocity error 


R-27940 



ERROR IN AV (fps) 


Figure 2.3-1 Error in a AV Pulse Count 

Due to Quantization 
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due to quantization is independent from sample to sample, the 
AV pulse count errors are correlated from sample to sample. 

Due to the nongaussian nature of this error source, a simpli- 
fied monte carlo simulation tool is used to support the evalua- 
tion of its contribution to calibration and alignment errors 
(Ref. 20). 


The quantization level for the resolvers is 20 sec. 
Because dithering rates are applied to all axes along which 
successive resolver measurements are made, this quantization 
error may be represented as an uncorrelated measurement noise 
with 


R n = (20) 2 /12 = 33.3 (sic) 2 


(2.3-5) 


The vehicle motions model was specified in detail in 
Ref. 1. The equivalent state space representation of this 
effect is 


- 


• 


- 




X., 


A 0 


X 


K 

r 

— w 

II 



-w 



N 

X 


B C 


K 


.0. 

rJ 

> 

1 




-V 




L J 


- 


- _ 

1 



where 


(2.3-6) 


*w 


North wind velocity 
West wind velocity 
West wind acceleration 


(2.3-7) 




North 1MU displacement 
North 1MU velocity 
West 1MU displacement 
West IMU velocity 
Vertical IMU displacement 
Vertical IMU velocity 


(2.3-8) 
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N is a white, gaussian random process and A, B, C and K are 
defined in Ref. 1. Since in the OFT mission both the Hangar 
and the Preflight calibration are performed on the launch pad, 
the same vehicle motion model is used for both cases. The I MU 
velocity induced by vehicle motions is modeled as a measurement 
error source for the LSF filter since the filter uses successive 
velocity measurements to derive 1MU accelerations and accelera- 
tion rates. 


2.3.3 Reliability of IMU Error Model and Data Base 

The majority of the error terms summarized by Tables 
2.3-2 through 2.3-4 are based on limited test data. Often the 
available test data may be from IMUs which do not incorporate 
the latest hardware and software modifications. This lack of 
sufficient representative test data may help to explain the 
fact that major changes of IMU error models have repeatedly 
been developed, often accompanied by extensive cal/align 
modifications. These modifications generally consist of 
costly and/or time consuming hardware and software fixes. A 
few examples are: 


• Two changes in the accelerometer scale 
factor asymmetry model, i.e., from the 
model used in the ALT software to that 
used in the OFT (Ref. 3) software to 
that recognized by Ref. 7 

• Belated recognition that accelerometer 
nonlinearities (Groups 18 and 19) are 10 
times larger in high gain in low gain 

• A recently proposed gyro randomness model 

• Recent identification of the accelero- 
meter double dead band problem 

• Several model and cal/align changes 
associated with gyro heading and 
accelerometer attitude sensitive effects. 
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These changes represent model improvements which are 
a natural consequence of the growing amount of applicable test 
data. The following questions, however, are worth asking: 


• How many additional model improvements 
are still to be discovered, and what 
would be their impact? 

• Is the "right kind" of test data being 
generated to help uncover potentially 
serious model (or IMU) deficiencies? 

• For any additional model deficiencies 
identified, is it really necessary to 
make hardware or software modifications? 


These questions are highly relevant because any assessment of 
IMU cal/align performance is ultimately limited by how well 
the assumed IMU error model matches the actual hardware error 
mechanisms . Recommendations for minimizing the impact of any 
IMU model deficiencies encountered subsequent to this cal/align 
analysis have been discussed with JSC. 


2. A OFT CALIBRATION AND ALIGNMENT PERFORMANCE 

The covariance analysis program for he IMU cal/align 
software includes all major changes to Hangar Coarse Alignment, 
Hangar Cal A, Preflight Coarse Alignment, Preflight Cal A, and 
Gyrocompassing as reflected by the 15 December 1978 release of 
the OFT Level C FSSR (Ref. 6). The error budget presented herein 
includes all error sources shown in Table 2.3-1 and discussed 
in Section 2.3 and Ref. 9. The results shown serve to verify 
the correctness of the cal/align equations specified in the 
FSSR and identify the dominant error mechanisms for all cali- 
brated IMU parameters. 
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The error budgets for Hangar Coarse Alignment, Hangar 
Cal A, Hangar Cal B, and Hangar Cal C are given in Section 

2.4.1, and for Preflight Coarse Alignment, Preflight Cal A, 
Gyrocompassing, and Velocity and Tilt Initialization in Section 

2.4.2. In addition, in Section 2.4.3, a more detailed error 
budget is given for the accelerometer attitude sensitive 
errors because these are dominant contributors throughout the 
entire calibration and alignment procedure. 


2.4.1 Hangar Calibration Performance 


The Hangar cal/align sequences have the following 
functions : 


• Hangar Coarse Alignment - Estimate the IMU 
case orientation. Table 2.4-1 displays 
both the rms level error and the azimuth 
error of this estimate. 

• Hangar Cal A - Calibrate a total of 35 
gyro and accelerometer error parameters. 
These parameters are organized into 10 
different groups. Table 2.4-1 displays 
the rms value of the calibration error 

of all parameters contained in each group. 
The corresponding column headings are 
identified by codes similar to those 
used in Table 2.3-1 (DF=gyro bias drifts, 
etc . ) . 

• Gangar Cal B - Calibrate 4 resolver 
error parameters and 2 misalignments of 
the compensated accelerometer frame 
with respect to the platform Z (azimuth 
rotation) axis. These parameters are 
organized into two groups. Table 2.4-1 
displays the rms value of the calibration 
error of all parameters contained in 
each group. The groups are accelerometer 
and gyro misalignments, and resolver 
offsets and gimbal nonorthogonalities, 
and the results are headed by 'GROUP 12' 
and 'GROUP 13', respectively. 
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CONTRIBUTING ERROR SOURCES 


CALIBRATED IMU ERRORS 

Gyro Bias Drift 

Gyro Hass Unbalance 

Gyro Torquer Scale Factor 

Gyro Misalignment 

Accel. Bias 

Accel. Scale Factor 

Accel. Bias 

Accel. Scale Factor 

Accel. Scale Factor Asymm. 


Low Gain 
Low Gain 
High Gain 
High Gain 
Correlated 
Uncorrelated 


11. Accel. Nonorthogonality 

12. Accel, and Gyro Misalignment 

13. Resolver Of fset , Gimbal Nonorthogonality 

II. NONCAL I BRATED IMU ERRORS 

14. Gyro Anisoelastici ty 

15. Gyro Output Axis Mass Unbalance 

16. Gyro Attitude Sensitive Drift 

17. Gyro Thermal Transient Drift 

18. Accel. Nonlinearity - 2nd Order 

19. Accel. Nonlinearity - 3rd Order 

20. Accel. Attitude Sensitive Bias 

21. Accel. Attitude Sensitive Scale Factor 

22. Outer Roll Offset and Misalignment 

23. Resolver Harmonic 

III. RANDOM AND QUANTIZATION ERRORS 


24. Gyro Randomness and Quantization | 

25. Accel. Randomness and Quantization i 

27. Vehicle Motion ) 

26. Resolver Randomness and Quantization 


MAGNITUDE 


0.035 sec /sec 
<0.300 sec/sec /g 
400 ppm 
120 Sefc 

<200 pg 
200 pg/g 
<200 pg 

200 pg/g 
200 Mg/fi 
30 Mg/g 
60 £eb 
60 sec 
<200 sec 


<0.025 sec/sec/g 
0.005 sec/sec/g 
<0.008 sec/sec 
0.01 sec/sec 
<150 Mg/gf 
< 50 ue/e J 


i 50 Mg/g 
<7.1 Mg 

< 22 Mg/g 
200 sec 

< 20 sec 


C*e Table 2.3-3 
See Table 2.3-3 


COARSE ALIGNMENT 


AZIMUTH 

sec 




II 


RSS OF ERROR GROUPS 2-27 


PERFORMANCE SPECIFICATION 


* These error budget contributions are negligible. 

* In low gain the error is 0 ppm, in high gain the error is 30 ppm. 

* 

Excluding contributions from Error Group 26 which are not significant. 

1. Determined from limits of tilt removal mechanization during calibration 

2. Value recommended by Singer-Kearfott . 

3. 30 sec for gimbal nonorthgonality . 

4. N.S. * Not Specified 



129.0 



4.9 

117.2 

3.8 

— 

28.9 

28.9 


23.7 

13.3 

13.3 

217.1 

784.1 

1600 1 

IBS 



< 3Too3 > 

0 003 


0.002 <Eii3> 
o.ooi 


< 3Too5 ) 14,. i <£ 


0.014, I 86.3 


0.025 200 6 


OI4D.OUT PRAMfi? 
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TABLE 2.4-1 

HANGAR CALIBRATION ERROR BUDGET 
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• Hangar Cal C - Determine the discrepancies 
between accelerometer and resolver indi- 
cation of platform level misalignments 
and store the difference of these discrep- 
ancies between each pair of IMU’s. These 
differences will be computed once again 
during Gyrocompassing and, when compared 
against the stored differences, constitute 
the Gyrocompass Goodness Test. 

Table 2.4-1 presents the performance results (both 
error budgets and overall performance) for these cal/align 
sequences. For reference, the magnitude of the largest error 

•L 

source within each error source group is given. The remaining 
numbers in Table 2.4-1 represent the root -sum- square (rss) of 
the contribution from all scalar error sources contained in 
each error source group, specified at the left of the corre- 
sponding row, to the calibration or alignment error of the 
error parameter group specified at the top of each corre- 
sponding column. As an example, the largest of the 11 gyro 

2 

anisoelasticities (Group 14) is 0.025 sec/sec/g . The rss 
effect of these 11 errors on gyro drift calibration (DF) is 
0.005 sec/sec per axis. 

For reasons of clarity, error contributions which are 
negligible are indicated by a dash (-) in Table 2.4-1. These 
values are all smaller than (most are much smaller than) 10% 
of the rss value of all error sources given in a particular 
column. Furthermore, the rss of all encircled numbers in each 
column is equal to at least 95% of the rss total for that 
column. Thus the encircled numbers represent dominant contri- 
butions t.o cal/align errors of a parameter group. A summary 
row near the bottom of Table 2.4-1 shows the combined rss 
effect of al 1 error groups, i.e., the projected cal/align 


Tables 2.3-2 to 2.3-4 give » magnitudes for all error sources. 
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performance. These results should be contrasted with the 
cal/align performance specification indicated in the last row 
of Table 2.4-1. 

Level 1MU case alignment errors at the end of Hangar 
Coarse Alignment are dominated, as expected, by resolver 
offset, gimbal nonorthogonality, and accelerometer and gyro 
misalignment. The azimuth coarse alignment error is dominated 
by west gyro bias drift and, to a lesser degree, by azimuth 
resolver offset, west gyro attitude sensitive drift (approxi- 
mately 0.008 sec/sec rms), and gyro misalignment toward West. 
Comparison of the last two rows in Table 2.4-1 indicates that 
the specified Hangar Coarse Aljjpunent accuracy appears to be 
achievable . 

The present Coarse Alignment results exhibit two 
significant deviations from those reported in Refs. 4 and 5 
for the algorithm defined in the December 1976 OFT Level C 
FSSR (Ref. 3). In the current Coarse Alignment mechanization, 
the IMU platform is earth rate compensated, and, therefore, 
gyro misalignment gives rise to a significant azimuth alignment 
error. In the earlier mechanization (Ref. 3) the platform was 
kept inertial, and, therefore, gyro misalignment had no signif- 
icant effect. Also, the effect, of gyro thermal transient 
drift is larger in Table 2.4-1 than in Ref. 4 because the 
magnitude assigned to this error source in the error model was 
increased (see Section 2.3.2). 

In Hangar Cal A , the dominant contributors to gyro 
bias and mass unbalance calibration error are noncalibrated 
gyro errors, primarily the attitude sensitive gyro drift. 

Gyro randomness and quantization also contribute significantly 
to mass unbalance calibration error, with the greatest impact 


2-26 


THE ANALYTIC SCIENCES CORPORATION 


Vr 

being on the input axis mass unbalance error . The gyro scale 
factor calibration error is dominated by second-order acceler- 
ometer scale factor nonlinearity. Accelerometer calibration 
errors are dominated by the effects of the uncalibrated attitude 
sensitivity of the accelerometer scale factor, discussed 
further in Section 2.4.3, and the second-order accelerometer 
nonlinearity. 

Comparison of the last two rows in Table 2.4-1 indicates 
that, except for low gain accelerometer bias and accelerometer 
scale factor asymmetry, the specified Hangar Cal A calibration 
accuracy appears to be achievable . The primary reason for 
noncompliance lies with accelerometer scale factor attitude 
sensitive errors. 


The present Hangar Cal A results exhibit several 
significant deviations from those reported in Ref. 4 for the 
cal/align mechanization described in the December 1976 OFT 
Level C FSSR. They are: 

• The effects of gyro attitude sensitive 
drift are reduced because the magnitude 
of this error source was decreased (see 
Section 2.3.2). 

• Gyro thermal transient drift is no 
longer a dominant error source; this 
fact owing to the reduction in the time 
constant of this error source (see 
Section 2.3.2). 

• Accelerometer second- and third-order 
nonlinearity has become a more important 
error source for those IMU calibration 
parameters which are estimated with 
accelerometers in high gain; this is due 
to the ten-fold increase of this error 


The contribution of gyro randomness and quantization to input 
axis mass unbalance error is 0.012 sec/sec/g. The smaller 
0.008 sec/sec/g is the rms contribution to all mass unbalance 
errors . 
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source in high gain (see Section 2.3.2). 
This effect is particularly noticeable 
with gyro torquer scale factor and high 
gain accelerometer bias calibration errors. 

• The dominant effect of accelerometer 

scale factor attitude sensitive errors 
on accelerometer calibration has increased 
even further, in spite of the fact that 
the magnitudes used for several of these 
error sources are smaller in the current 
OFT cal/align performance analysis. 

This result is not surprising, considering 
that the former mechanization was specially 
designed to be insensitive to this error 
source, whereas the current OFT cal/align 
mechanization is not. More detailed 
results on this error source are given 
in Section 2.4.3. 


In Hangar Cal B , resolver randomness errors and quan- 
tization contribute significantly to all calibration errors, 
although resolver offset and gimbal nonorthogonality calibra- 
tion errors are dominated by uncalibrated resolver harmonic 
errors. The results in Table 2.4-1 do not differ significantly 
from those obt .ined in Refs. 4 and 5. Comparison of the last 
two rows in Table 2.4-1 indicates that the specified Hangar 
Cal B calibration accuracy appears to be achievable . 

Hangar Cal C errors are dominated by resolver har- 
monics and, to a lesser degree, by accelerometer second-order 
nonlinearity and attitude sensitive scale factor errors. The 
results in Table 2.4-1 do not differ significantly from those 
obtained in Refs. 4 and 5. There is no performance specifica- 
tion for Hangar Cal C. The Hangar Cal C results are interme- 
diate and will be used for the goodness test in gyrocompassi ng. 
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In summary, much of the Hangar Calibration error buck- 
et appears to be dominated by the uncalibrated accelerometer 
scale factor attitude sensitivity . In view of the importance 
of this error source, further discussions of its impact are 
provided in Section 3.4.3. The gyro scale factor and high gain 
acceleromter bias calibration errors are dominated by second- 
order accelerometer nonlinearity . The gyro bias and mass unbal- 
ance calibration errros are dominated by gyro attitude sensitive 
drift and, in the case of input axis mass unbalances, also by 
gyro randomness. 


2.4.2 Preflight Calibration and Alignment Performance 


The Preflight cal/align sequences have the following 
functions : 


• Preflight Coarse Alignment - Estimate 
the 1MU case orientation. Table 2.4-2 
displays both the rms level error and 
the azimuth error of this estimate. 

• Preflight Cal A - Calibrate a total of 

1& gyro and acceleromter error parameters. 
Table 2.4-2 displays the rms values of 
all error groups calibrated in this se- 
quence. The other error groups which 
were calibrated in Hangar Cal A but not 
in Preflight Cal A are also included for 
reasons of clarity. 

• Gyrocompassing - Estimate the platform 
azimuth misalignment and perform a good- 
ness test. Table 2.4-3 displays the rms 
values of the level and azimuth errors 
expected during the goodness test, as 
well as the rms value of the azimuth 
alignment error. 

• Velocity and Tilt Initialization - Perform 
a final level alignment of the IMU and 
provide a transition period before lift- 
off. Table 2.4-3 exhibits the level and 
azimuth misalignment error at the end of 
Velocity and Tilt Initialization. 
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CONTRIBUTING ERROR SOURCES 
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TABLE 2.4-2 

PREFLIGHT CALIBRATION ERROR BUDGET 



F°LDOUX FRAME Z. 
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TABLE 2.4-3 

PREFLIGHT FINE ALIGNMENT ERROR BUDGET 
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Tables 2.4-2 and 2.4-3 summarize Che performance 
results (both error budgets and overall performance) for these 
cal/align sequences. Most entries in these tables are analogous 
to those of Table 2.4-1. In addition. Table 2.4-2 shows the 
largest turn-off turn-on increment magnitude in each error 
source group. These error increments are added in an rss 
sense to the residual calibration errors at the end of Hangar 
calibration. The salient features of the Preflight cal/align 
performance evaluation results are discussed below. 

Level IMU case alignment errors in Preflight Coarse 
Alignment are dominated primarily by second-order accelerometer 
nonlinearity and by resolver harmonic, offset, and nonorthogo- 
nality: while azimuth aligment errors are dominated by the 
turn-off turn-on increment in gyro bias drift and by gyro 
attitude sensitive drift. Note that the contribution of 
accelerometer second-order nonlinearity, accelerometer scale 
factor attitude sensitivity, and gyro attitude sensitive drift 
to Preflight Coarse Alignment errors is greater than the 
corresponding contribution in Hangar Coarse Alignment, due to 
their influence on Hangar Cal A calibration errors (Table 
2.4-1). Comparison of the last two rows of Table 2.4-2 
indicates that the specified Preflight Coarse Alignment 
accuracy appears to be achievable . 

The present Coarse Alignment results exhibit two 
significant deviations from those reported in Refs. 4 and 5 
for the algorithm defined in the December 1976 OFT Level C 
FSSR (Ref. 3). The level alignment accuracy is now dominated 
by the second-order accelerometer nonlinearity, owing to the 
ten- fold increase of this error source. The azimuth alignment 
error increased substantially, owing to the increase in the 
gyro bias drift turn-off turn-on instability. 


Tables 2.3-2 to 2.3-4 give the magnitudes for all error sources. 
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The Preflight Cal A error budget is similar to the 
Hangar Cal A error budget, with two noticeable differences: 

First , the IMU errors which are calibrated in Hangar Cal A but 
not recalibrated in Preflight Cal A (f in Table 3.4-2) show 
degradation due to turn-off turn-on instabilities. In most 
such cases these instabilities dominate the error budget. 

Second, since accelerometer scale factor asymmetry is not 
recalibrated, its turn-off turn-on instability contributes 
directly to accelerometer bias calibration error. 

Comparison of the last two rows in Table 2.4-2 indi- 
cates that except for low and high gain accelerometer bias, 
the specified Preflight Cal A accuracy appears to be achievable . 
The primary reason for noncompliance lies with the accelerometer 
scale factor attitude sensitive error. 

The IMU remains on after Preflight Cal A, which insures 
against parameter shifts due to power turn-off turn-on. The 
salient features of the Preflight Alignment performance evalua- 
tion results are discussed below. 

The Gyrocompass Goodness Test error budget is pre- 
sented in Table 2.4-3 in the form of rms relative level and 
azimuth IMU misalignment errors. The level channels of the 
goodness test are affected primarily by turn-off turn-on 
instabilities. In the azimuth channel (with the baseline 
assumption that perfect compensation is achieved for gyro atti- 
tude sensitive drifts, resolver offset, and gimbal nonorthogon- 
ality), second-order accelerometer nonlinearity and gyro ran- 
domness would be the dominant error sources. However, the 
accuracy of the azimuth axis goodness test is extremely sensi- 
tive to gyro attitude sensitive drift compensation errors . For 
example, if no compensation were applied, the rms accuracy of 
the relative azimuth misalignment used in the goodness test 
would degrade to approximately 440 sec. Comparison of the 
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last two rows in Table 2.4-3 indicates that, assuming perfect 
gyro attitude sensitive drift compensation, the residual values 
of the goodness test are within the threshold limits specified 
for that test . 

The Gyrocompass azimuth alignment error budget is 
similar to the azimuth channel of the goodness test, with two 
notable differences: First , resolver offsets and gimbal non- 
orthogonalities do not play a role since resolver measurements 
are not used for azimuth alignment. Second , the statistical 
effects of most other error sources are reduced by a factor of 
■[l since the goodness test is a comparison of two IMU's while 
azimuth alignment utilizes measurements from only one 1MU. 
Randomness and quantization errors are significant error con- 
tributors in azimuth alignment. As in the azimuth goodness 
test, there is a strong sensitivity to errors in the gyro 
heading sensitive drift compensation , yielding a degradation 
in azimuth alignment error of about 313 sec in the case of no 
compensation. In order to achieve the specified azimuth align- 
ment accuracy of 60 sec, approximately 85% of the attitude 
sensitive drift error at the gyrocompass positions must be 
removed via special calibration. To insure satisfactory system 
performance, this calibration should be checked at regular 
intervals . 


Comparison of the last two rows of Table 2.4-3 indi- 
cates that the performance specification for azimuth alignment 
appears to be achievable. This conclusion, however, assumes 
that near perfect compensation for gyro heading sensitive drifts 
can be achieved. 


The Velocity and Tilt Initialization error budget is 
also presented in the form of rms level and azimuth alignment 
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errors. Level alignment errors seem to be dominated by the 
turn-off turn-on instability in accelerometer misalignments 
and by accelerometer attitude sensitve errors. The azimuth 
alignment error budget is similar to that in Gyrocompassing, 
and the same comments apply with regard to gyro heading sensi- 
tive drift compensation. Comparison of the last two lines of 
Table 2.4-3 indicates that the performance specification for 
final level alignment appears achievable 

In summary , at the end of preflight calibration and 
alignment, a review of calibration and alignment errors aris- 
ing from all error sources indicates that: 


• The specified calibration accuracies for 
low and high gain accelerometer biases 
and for accelerometer scale factor asym- 
metry are not met on the basis of the 
current 1MU error model. The primary 
reason for noncompliance lies with the 
accelerometer scale factor attitude sen- 
sitive errors. 

• The performance specification for all 
other 1MU cal/align error appears to be 
acheivable . 

• Gyro calibration errors are dominated by 
turn-off turn-on instabilities. 

• Accelerometer calibration errors are 
dominated by attitude sensitive scale 
factor errors (discussed further in Sec- 
tion 2.4.3), second-order accelerometer 
nonlinearity, and turn-off turn-on insta- 
bility. 

• The primary error sources in the baseline 
Preflight Alignment evaluation are: 
resolver instabilities in the goodness 
test, gyro randomness and quantization 

in the azimuth measurement, and acceler- 
ometer misalignment instability and non- 
linearity in the level axis measurements. 
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• The azimuth alignment measurements are 

extremely sensitive to errors in attitude 
sensitive gyro drift compensation, 
yi eld ing rms azimuth alignment errors of 
5 min in each IMU or 7 mxn between two 
IMUs if no compensation is applied. 
Approximately 85% of this error must be 
removed via periodic calibration in order 
to achieve the azimuth alignment specif i- 
cation. 


2.4.3 Effect of Accelerometer Scale Factor Attitude 
Sensitivity 


Due to the importance of the accelerometer scale 
factor attitude sensitivity (error Group 21) in the perform- 
ance and evaluation results presented in Sections 2.4.1 and 
2.4.2, a more detailed error budget for this error group is 
presented here. 

A detailed error budget, showing the contribution of 
each accelerometer scale factor attitude sensitivity parameter 
to each accelerometer calibration error in Hangar Cal A is 
shown in Table 2.4-4. Each row of the table represents an 
individual error in Group 21 , and each column represents an 
individual accelerometer calibration error in one of Groups 6 
to 10. (The high gain results are identical to the low gain 
results. ) 


As can be seen in Table 2.4-4, the accelerometer bias 
and asymmetry calibration errors are affected to much greater 
degree than are the scale factor calibration errors. The 
largest error contributions come from the X accelerometer roll 
gimbal sensitivity, followed by the Z accelerometer azimuth 
gimbal sensitivity. 
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TABLE 2.4-4 

DETAILED ERROR BUDGET FOR 
ERROR GROUP 21 IN HANGAR CAL A 
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Because of the overwhelming effect of the accelerometer 
attitude sensitive errors it was decided (see Ref. 7) to dynam- 
ically compensate these errors within the 1MU software. Given 
this situation, it is of interest to know how small the residual 
attitude sensitivity errors (after compensation) have to be in 
order that this error source ceases to be a dominant contributor 
to cal/align errors. To help answer this question, the error 
budget contributions in Table 2.4-4 were normalized by dividing 
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these by the magnitude of the corresponding error sources; the 
latter are also shown in Table 2.4-4. The result of this oper- 
ation is shown in Table 2.4-5 which contains the factors by 
which an error source is amplified to determine its contribution 
to a calibration error. A similar table was constructed for 
the accelerometer bias attitude sensitivity errors which, as 
per Ref. 7, will a 1 so be compensated within the IMU software. 
Table 2.4-6 contains .he corresponding error amplification 
factors. 


In Table 2.^-5 it can be seen that error amplification 
factors as high as 4.4 and 4.7 occur several times. The highest 
factor in Table 2.4-6 is 3.4. Information on these amplifica- 
tion factors was provided to RI to help them specify accuracy 
requirements for the dynamic compensation of accelerometer 
attitude sensitivity. 

The amplification factors shown in Table 2.4-5 and 
2.4-6 depend strongly on the particular platform positions 
used for calibration purposes. The fact that many of these 
factors are so large is an indication that the current cal/ 
align mechanization is highly sensitive to accelerometer 
attitude sensitive errors . 

2.4.4 Needs for IMU Error Model Improvement 

The results discussed in Section 2.4 clearly identif 3 » 
a number of error sources which are important contributors to 
cal/align performance. The mathematical models used for these 
error sources are described in Section 2.3. However, as 
pointed out in Section 2.3.3, the majority of these models are 
based on limited test data and, therefore, may not accurately 
reflect the actual size and characteristics of the corresponding 
error mechanisms. Thus, future projections of cal/align perform- 
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TABLE 2.4-5 

ERROR AMPLIFICATION FACTORS FOR 
ERROR GROUP 21 IN HANGAR CAL A 


T-3274 



ance can be improved primarily by seeking better models for those 
error sources which present study has found to be most limiting . 
The latter includes the following: 


• The residual error (after dynamic compensa- 
tion) of accelerometer scale factor attitude 
sensitivity (Group 21) , because of its ~ 
dominating effect on acceleromter calibra- 
tion errors 
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TABLE 2.4-6 

ERROR AMPLIFICATION FACTORS FOR 
ERROR GROUP 20 IN HANGAR CAL A 
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• Second-order accelerometer nonlinearity 
TGroup ll) , both in high and low gain, 
because of its dominati ' effect on 
accelerometer calibration errors and 
gyrocompass level and azimuth alignment 

• Turn-off turn-on errors of input axis 
gyro mass unbalance (Grou^ 2) , because 

of their dominating effect on gyrocompass 
azimuth alignment 
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• Turn-off turn-on instabilities of 
accelerometer and gyro misalignments 
{Group 12) , because of their qegradation 
of gyrocompass level alignments 

• Gyro randomness and quantization (Group 24) , 
because of their degradation of gyro 

input axis mass unbalance calibration 
and gyrocompass azimuth alignment 

• The accuracy and stability of the gyro 
heading sensitive drift compensation , 
because of its impact on gyrocompass 
azimuth alignment. 


It is recommended that the following steps be taken 
in relation to the above mentioned dominant error sources: 


• Generation and evaluation of test data 

to verify and/or improve the corresponding 
error models 

• Evaluation of the impact of the improved 
error models on cal/align performance 

• Evaluation of the effect of dominant 
error sources on orbit insertion accuracy. 


2.5 SUMMARY AND CONCLUSIONS 

A detailed performance analysis of. the Space Shuttle 
alibration and alignment software has been completed. The 
software evaluated corresponds to the 15 December 3978 issue 
of the OFT Level C FSSR. The IMU error models used in this 
study have been updated with the most recent information avail- 
able on IMU performance. 

2.5.1 Summary of Findings 

The performance evaluation of the cal/align software 
has included two efforts. They are: 
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• Review of the error models used to repre- 
sent all potentially important IMU error 
mechanisms 

• Development of an error budget showing 
the contribution of all modeled error 
sources to calibration and alignment 
errors . 


The principal findings of both efforts are summarized below. 


Analytical evaluation of IMU cal/align performance is 
limited by the completeness and accuracy of the models used to 
represent IMU errors. For this reason the present study has 
included a complete review of IMU error models. The major 
findings of this effort are as follows: 


• The error models for accelerometer scale 
factor asymmetry, accelerometer nonlineari- 
ties, and gyro randomness, as well as 
several error model parameters, were modi- 
fxed to reflect information obtained 

from recent IMU test data. The final 
error model used by TASC in this study 
was reviewed by SKD personnel. 

• Many of the above error model improve- 
ments have a significant effect on cal/ 
align performance. Since the need for 
these model improvements has been identi- 
fied only recently, there is some concern 
that there may be other important model 
deficiencies which have not yet been 
discovered. 


Recommendations for further improvements of IMU error models 
and, hence, for further improving the confidence of the cal/ 
align performance projections are given in Section 3.2. 

The evaluation of the calibration and alignment 
performance was based upon a covariance analysis of all 
currently known IMU errors. The major findings of this effort 
are as follows: 


2-42 


THE ANALYTIC SCIENCES CORPORATION 


• The specified calibration accuracies for 
low and high gain accelerometer biases 
and for accelerometer scale factor asym- 
metry are not met. The primary reason 
for noncompliance lies with the acceler- 
ometer scale factor attitude sensitive 
errors . 

• The performance specification for all 
other 1MU cal/align errors appears to be 
achievable. 

• Gyrocompass azimuth alignment is extremely 
sensitive to errors in attitude sensitive 
gyro drift compensation, yield in g rms 
azimuth alignment errors of 5 fnih in 

each 1MU or 7 min between two IMU's if 
no compensation is applied. Approximately 
85% of this error must be removed via 
calibration in order to achieve the azimuth 
alignment specification. 

• Gyro calibration errors are dominated by 
turn-off turn-on instabilities. 

• Accelerometer calibration errors are 
dominated by attitude sensitive scale 
factor errors, second-order accelerometer 
nonlinearity, and turn-off turn-on insta- 
bility. 

• The primary error sources in the base- 
line Preflight Alignment evaluation are: 
resolver instabilities in the goodness 
test, gyro randomness and quantization 
in the azimuth measurement, and acceler- 
ometer misalignment instability and 
nonlinearity in the level axis measure- 
ments . 


The accelerometer attitude sensitive scale factor 
errors will be dynamically compensated within the 1MU software 
(c.f. Ref. 7). Test data on the compensated IMU's should be 
generated and evaluated to develop an error model for the resi- 
dual attitude sensitive errors. Based on the detailed cal/align 
error budget generated for the study reported here, TASC has 
provided NASA and Rockwell International with information 
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necessary to specify the minimum required performance of the 
above mentioned compensation. 

2.5.2 Recommendations 

The majority of the Space Shuttle 1MU error models 
are based on limited test data and, therefore, may not accu- 
rately reflect the actual size and characteristics of the 
corresponding error mechanisms. Thus, estimation of cal/ 
align performance under mission conditions can be improved 
primarily by seeking better models for those error sources 
which the present study has found to be most limiting . The 
latter include: 


• Accelerometer scale factor attitude 
sensitivity after dynamic compensation 

• Second-order accelerometer nonlinearity 

• Turn-off turn-on instabilities of input 
axis gyro mass unbalance and accelerometer 
and gyro misalignments 

• Gyro randomness and quantization 

• Gyro attitude sensitive drift compensa- 
tion for the gyrocompass positions. 


It 

in relation 


is recommended that the following steps be taken 
to these dominant error sources: 


• Generation and evaluation of test data 
to verify and/or improve the corre- 
sponding error models 

• Evaluation of the impact of the improved 
error models on cal/align performance 

• Evaluation of the effect of dominant 
error sources on orbit insertion 
accuracy. 
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3. POSTFLIGHT IMU ERROR RECOVERY FOR APPROACH AND 

LANDING MISSION PHASES 


Prior to the initial Space Shuttle Approach and Landing 
Tests (ALT), NASA's experience with the Shuttle IMUs was gathered 
from bench tests and hardware/ software integration tests. The 
alignment phase of the bench tests is designed to evaluate IMU 
alignment under laboratory conditions in absolute terms; however, 
the calibration phase is designed only to evaluate the repeat* 
ability of the calibrated IMU parameters. Atmospheric flight 
tests such as ALT provide an opportunity to evaluate both the 
navigation performance achievable with the IMU and the true 
calibration and alignment uncertainties in a flight environment. 
This chapter provides assessment of the potential accuracy of 
postflight IMU error recovery for atmospheric flight tests 
similar to ALT. 

There are two likely benefits of postflight IMU error 
recovery. The most obvious is to validate the IMU performance 
(both hardware and software). A second benefit of equal impor- 
tance is that postflight recovery of inflight IMU alignment 
errors would provide an accurate attitude reference for evalu- 
ating (and calibratiig) other Shuttle subsystems such as the 
air data system. Postflight filtering of flight test data is 
sufficient for the IMU performance evaluation; however, smooth- 
ing of alignment errors is desirable for generating the attitude 
reference. Although the analysis presented in this chapter is 
restricted to an evaluation of postflight filtering accuracy, 
it is possible to make strong inferences on smoother performance 
based on these results. The ALT mission scenario and system 
error models used in this study are discussed in Section 3.1. 

The performance results are presented in Section 3.2. 
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3.1 MISSION SCENARIO AND SYSTEM ERROR MODELS 

There are two principal ingredients of a successful 
performance test program: 

• The test design must be adequate to excite 
the important error mechanisms. 

• The test instrumentation must be accurate 
enough to measure the resulting perform- 
ance errors. 

Because the dynamic environment for the approach and landing 
mission phases is benign relative to the environment for the 
ascent and reentry phases, some 1MU error sources (in particular, 
g-sensitive IMU errors) are not observable from ALT test data. 
Nonetheless, it is possible to select an ALT trajectory to 
maximize post flight IMU error recovery. In particular, the 
duration of the test flight should be sufficient for gyro drifts 
and acceleromter biases to be estimated, and the trajectory 
should have enough turning maneuvers to make azimuth misalign- 
ments observable. 

The ALT test scenario selected for the study satisfies 
both of the above criteria. The ground track for the trajectory 
is presented in Fig. 3.1-1, and a mission time line is given 
in Table 3.1-1. The trajectory closely corresponds in both 
geometry and time-line to the two-abort captive trajectories 
described in Ref. 21. The principal differences are that, for 
simplicity, the takeoff and landing were not simulated (dashed 
curves in Fig. 3.1-1), and the remainder of the trajectory is 
at constant velocity (500 fps) and altitude. Measurements taken 
during takeoff and landing would aid error recovery somewhat, 
but would not significantly alter the performance results pre- 
sented in Section 3.2. 
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Figure 3 


••tsoii 



.1-1 Ground Track for Sample ALT Trajectory 
TABLE 3.1-1 

TIME LINE FOR SAMPLE ALT TRAJECTORY 


T-0865 


EVENT 

TIME 

Preflight 


Cal /Align 

-60 min 

Turn 1 Exit 

0 

Turn 2 Enter 

$ min 

Exit 

9 min 50 sec 

Turn 3 Enter 

29 min 10 sec 

Exit 

33 min 40 sec 

Turn 4 Enter 

48 min 40 sec 

Exit 

53 min 10 sec 

Turn 5 Enter 

68 min 10 sec 

Exit 

72 min 40 sec 

Turn 6 Enter 

88 min 40 sec 

Exit 

91 min 10 sec 

Turn 7 Enter 

98 min 40 sec 


3-3 



THE ANALYTIC SCIENCES CORPORATION 


For study purposes, a generalized test instrumentation 
model that corresponds roughly to photo theodolites in terms 
of measurement type and accuracy was selected. A position fix 
of approximately 30 ft accuracy per axis (one sigma) was assumed 
to be available at 30 sec intervals along the entire trajectory. 
A summary of the instrumentation error model is presented in 
Table 3.1-2. 

As previously mentioned, ALT trajectories are not 
sufficiently active to make g-sensitive 1MU errors (gyro mass 
unbalances, accelerometer scale factors, accelerometer asym- 
metries, etc.) observable. Hence, there is no significant 
advantage to including these error sources in either a filter 
or a truth model for ALT postflight error recovery. An excep- 
tion is the input axis mass unbalance for the z-gyro (vertical 
input axis). This error source is continually excited by grav- 
ity and is indistinguishable from the z-gyro bias drift for 
ALT. 


The IMU error model for the study is presented in 
Table 3.1-3. The only error sources included are those which 
are expected to be observable, except that the acceleromter 
scale factor errors were included to insure that the azimuth 
misalignment estimate did not become overly optimistic. 


TABLE 3.1-2 

ASSUMED INSTRUMENTATION ERROR MODEL FOR ALT 


MEASUREMENT TYPE 

m 

MEASUREMENT 

NOISE 

(la) 

Position 

15 ft 

25 ft 
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TABLE 3.1-3 

DATA BASE FOR IMU-RELATED ERROR SOURCES 


ERROR SOURCE 

STANDARD DEVIATION 

Initial Misalignments 


Azimuth 

80 sefc* 

Tilt 

20 sefc* 

Accelerometers 


Biases 

50 pg 

Scale Factors 

100 ppm 

Quantization 

1 cm/sec 

Gyros 


Bias Drifts 

0.015 deg/hr 

Mass Unbalances 

0.025 deg/ hr /g 


♦One hour prior to launch. 


The model reflects the nominal prelaunch calibration and align- 
ment performance (see Chapter 2) except that the acceleromter 
scale factor uncertainty has been increased to account for 
other effects such as asymmetries and nonlinearities. The 1MU 
errors at completion of the calibration and alignment sequence 
are modeled as uncorrelated except for a 0.7 correlation co- 
efficient between tilt misalignments and accelerometer biases. 

A typical ALT mission plan requires the 1MU calibra- 
tion and alignment sequence to be completed approximately one 
hour prior to takeoff. The effect of this delay is to intro- 
duce an increase in the 1MU tilt and azimuth misalignments 
due to the gyro drifts along the appropriate axes. The delay 
was accounted for in the simulation by increasing the tilt 
and azimuth misalignments to 40 sec and 120 sec, respectively, 
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and by adding the resulting correlations between alignment 
and gyro errors. 


3.2 ASSESSMENT OF POTENTIAL ERROR RECOVERY PERFORMANCE 

The results of a covariance analysis of an optimal 
postflight processor for the simulated ALT mission are pre- 
sented in Table 3.2-1 for three times: after the calibration 

and alignment sequence, at the start of the flight, and at the 
end of the flight. Because gyro random drift was omitted from 
the error model, the filtered gyro errors at the start of the 
flight are the same as at the end of the calibration and align- 
ment sequence. The accelerometer uncertainties are also con- 
stant over this interval, but this is consistent with the com- 
plete IMU error model in Chapter 2. 


Table 3.2-1 indicates that the ALT measurements are 
sufficiently accurate to yield a significant reduction in all 


TABLE 3.2-1 

STANDARD DEVIATIONS OF FILTERED IMU ERRORS 
FOR SAMPLE ALT MISSION 


T-0478 




ERROR SOURCE (la) j 

EVENT 

TIME 

(min) 

TILT 

MISALIGNMENT 

(dec) 

AZIMUTH 

MISALIGNMENT 

(Sec) 

lifflij 

VERTICAL AXIS 
GYRO DRIFT* 
(deg/lir) 

ACCELEROMETER 

BIAS 

(Ug) 

ACCELEROMETER 
SCALE FACTOR 
(ppm) 

Preflight 
cal /align 

-60 

20 

80 

0.015 

0.030 

50 

100 

Exit 
Turn 1 

0 

40 

120 

0.015 

0.030 

50 

100 

Enter 
Turn 7 

100 

11 

70 

0.002 

0.013 

40 

95 


•Includes input axis mass unbalance. 


*Gyro random drift is modeled as a random walk which would amount 
to only 0.006 deg/hr at the end of one hour (see Chapter 2). 

This would rss to increase the gyro bias drift from 0.015 deg/hr 
to 0.016 deg/hr. 
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the level axis errors (tilts, drifts for gyros in the level 
plane, and acceleromter biases) except for accelerometer 
scale factor errors. The improvement in the azimuth error 
estimates (azimuth misalignment and vertical gyro drift) is 
less dramatic, but the final uncertainty for these parameters 
is still smaller than it was immediately following the cali- 
bration and alignment sequence. 

Figure 3.2-1 provides a graph of the IMU misalign- 
ment errors as a function of time. As expected, the azimuth 
uncertainty is reduced by the horizontal accelerations accom- 
panying each turn. Whereas the tilt uncertainty is reduced 
to a steady state value almost immediately, the azimuth error 
requires 50 min of flight (three turns) to approach a steady 
state value. A closer inspection of the data suggests that a 
standard deviation of 60-70 sec may be the smallest azimuth 
error attainable with ALT measurements, primarily because of 
velocity errors attributable to accelerometer quantization. 



Figure 3.2-1 Filtered IMU Misalignments 

for Simulated ALT Mission 
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In addition to the fact that some errors (such as 
accelerometer scale factor) are unobservable for ALT because 
they generate small position errors in the ALT environment, 
other errors are nonseparable because two or more error 
sources have similar dynamic effects in a one-g environment. 
Gyro bias drift and input axis mass unbalance for the vertical 
axis gyro provide one example - the ALT instrumentation meas- 
ures the sum of the navigation error due to these two error 
sources and thus the error estimates generated from postflight 
analysis will be highly correlated. Similarly, the final tilt 
and accelerometer bias estimates will be almost perfectly cor- 
related. For the error models assumed in this study, the 
nonseparability of individual error sources is more of a factor 
in limiting postflight error recovery than the accuracy of the 
ALT instrumentation set . 

Because the gyro and accelerometer errors are con- 
stant, the filter performance uncertainties predicted in Table 
3.2-1 are the same as the smoother performance would be. Thus, 
the only estimates which would be improved by smoothing are 
those for the inflight misalignments - and, for reasons dis- 
cussed in the preceding paragraphs, it is doubtful that accura- 
cies much better than 60 sec and 10 sec can be attained for 
azimuth and tilt, respectively. Nonetheless, it appears that 
postflight processing of ALT measurements could maintain IMU 
alignment errors at a level smaller than the preflight misalign 
ment values . Even with the additional attitude uncertainty 
associated with the gimbal resolvers (on the order of 30 sec, 
see Chapter 2), this should provide a useable reference for 
estimating such attitude-dependent variables as Orbiter angle- 
of-attack and sideslip angle, and hence for calibrating the 
air data system. 
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4. ON ORBIT CALIBRATION OF SHUTTLE IMU 

ACCELEROMETER AND GYROSCOPE BIASES 


4.1 INTRODUCTION AND SUMMARY 

A simple calibration of the Shuttle IMU accelerometer 
and gyroscope biases is planned for orbital operations (Ref. 
3). In the very low sensed acceleration environment on orbit, 
only the instrument biases (as opposed to scale factor errors, 
etc.) are observable. The on-orbit calibration is an attempt 
to reduce the likely errors in compensating for these biases 
to values equivalent to, or better than, those of the last 
ground calibration, thus mitigating the effects of on-orbit 
bias instability. Reasonable target values for the accuracy 
of on-orbit calibration are thus 0.15 deg/hr (lo) gyro bias 
and 50 pg (lo) accelerometer bias (Ref. 22). This chapter 
examines the likely accuracy of the on-orbit calibration pro- 
cess . 


The calibrations for accelerometer biases and for 
gyro biases are quite different. The accelerometer calibra - 
tion takes only 320 sec, so that long-term instrument stabil - 
ity is of little importance to the calibration process itself . 
Also, accelerometer outputs are observed directly by accumula- 
ting net delta-velocity pulses, so that instrument measurement 
errors, except for quantization effects, are negligible. The 
accelerometer outputs are, however, very sensitive to dynamical 
disturbances (centripetal acceleration , venting, etc), and 
these disturbances impose lower limits on the accuracy of accel- 
erometer bias calibration. On the other hand, the duration of 
the gyroscope calibration is yet to be determined, but it is 
anticipated to be several hours. The gyro drift cannot be 
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expected to remain constant during this interval, because of 
instrument instabilities. While gyro outputs are not appreci- 
ably affected by dynamical disturbances, neither is gyro drift 
directly observable. A star-tracker-to-IMU-cluster alignment 
procedure must be executed at the beginning and end of the 
gyro calibration period to assess net gyro drift during the 
period, and alignment measurement errors are significant. Thus , 
gyro calibration is largely a tradeoff of instrument stability 
and measurement error; selecting a calibration period to balance 
the two will produce the best accuracy . 

An important consideration in gyro instrument sta- 
bility is a type of instability known as heading sensitivity. 
Heading sensitivity refers to variations in gyro drift rate 
with a change in cluster-to-IMU-case attitude (see Chafer 2). 
Also included under the heading sensitivity category are varia- 
tions in gyi drift rate with the directions the IMU gimbals 
are rotating and with the direction in which the gyro was last 
torqued for slewing. These effects have been observed in a 
one-g ground environment and it is considered likely that they 
will also exist in a zero-g environment. Accelerometers are 
also subject to some heading sensitive effects, but these effects 
are largely scale factor instabilities that are unimportant in 
the orbital environment. 

The following two sections (Sections 4.2 and 4.3) on 
accelerometer calibration and on gyro calibration consider 
these processes in more detail. The principal conclusions of 
this Chapter are presented in Section 4.4. 

• Both of the bias calibration processes 

appear feasible in their proposed forms; 
no more sophisticated or complex processes 
are required. 
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• With mild controls on dynamical disturbances: 

- Thrust to 2 lb. rms or less 

- Vehicle rotation rate to 0.15 deg/sec 
or less 

- Unpaired impulsive events to 20,000 
lb-ft/sec or less 

accelerometer bias calibration accuracy 
should be better than 20 |jg(lo). 

• The amount and effect of heading sensitive 
gyro drift is very difficult to predict 

on orbit. Without it, a 16 hr gyro calibra- 
tion period is optimal, and would be 
accurate to 0.007 deg/hr (lo). Any period 
from 7 hr to 36 hr would yield bias drift 
rate calibration accuracy of better than 
0.01 deg/hr (lo). With heading sensitive 
gyro drift, the accuracy of the bias 
calibration could degrade to 0.02 deg/hr 
(lo) or worse. 


4 . 2 ACCELEROMETER CALIBRATION 


In accelerometer calibration on orbit, net compen- 
sated accelerometer pulses accumulated over a time interval 
are divided by the interval length to estimate uncompensated 
accelerometer bias (Ref. 3). The procedure assumes that there 
is a negligible input to the accelerometers on orbit, so that 
all of the compensated output represents uncompensated bias 
effects. Errors in the calibration result from the quantiza- 
tion of the accelerometer output into digital pulses, from 
disturbing dynamical effects, and (to a lesser extent) from 
accelerometer attitude sensitivity. 

In the low gain mode used for the calibration, one 
accelerometer pulse represents one milli-g-sec of velocity 
change. An accelerometer calibration interval of 320 sec 
is contemplated, so that a one-pulse error would produce a 
bias estimation error of 3.125 pg. If the precise beginning 
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and ending tiroes of the calibration were selected without re- 
gard to the accelerometer pulse stream itself, and if the 
accelerometer produced pulses uniformly on orbit, then the rms 
bias estimation error due to quantization would be only about 
40% of the one-pulse error, or about 1.28 pg. Accelerometers, 
however, do not all behave ideally. Some instruments output 
pulses in groups of two or more with the same polarity, even 
for vanishingly low input accelerations. Given that the tar- 
get accuracy of accelerometer bias compensation is in the 50 
pg range, even moderately abnormal quantization effects will 
not likely prove to be a limiting factor. 

Disturbing effects due to vehicle motion are another 
matter. Any large disturbances would be detected by the un- 
reasonableness of the calibration test results themselves, and 
moderate disturbing effects localized in time to periods much 
shorter than 320 sec could be caught by a well designed test 
that used sub- interval data for reasonableness checks. Some 
disturbances will have to be controlled or monitored by other 
means . 


A list of likely disturbance . sources is: 


• Thrust due to gasses venting from the 
orbiter, possibly for water evaporation 
(for cooling) or for vehicle attitude 
control 

• Vehicle rotation rates, possibly due to 
star acquisition for concurrent gyro 
calibration 

• Impulsive events, possibly due to pay- 
load shifting or hatch movement. 


RMS venting thrusts are currently estimated to be in 
the 1 lb region (Ref. 23). Assuming a 200,000 lb Orbiter, 1 lb 
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of thrust gives 5 pg of acceleration, which translates direct- 
ly into a 5 |jg acceleromter bias error during calibration. 

This small value would be acceptable. Some exceptional venting 
thrusts can produce unacceptable errors. The use of main atti- 
tude control jets for major attitude changes car. result in a 
30 lb burst of thrust over a 2 min span. Such an event must 
not be allowed during accelerometer calibration, for it would 
produce a 56 pg error. In general, 1 lb-min of thrust-times- 
duration in accelerometer calibration will produce 0.94 pg of 
bias error. 

Vehicle rotations can produce sizeable acceleration 
errors because the IMU is located about 50 ft from the Orbiter 
center of mass. Centripetal acceleration is given by the famil- 
iar 


a = w 2 r , (4.2-1) 

where u is the rotation rate, and r is the perpendicular dis- 
tance from the rotation axis. Assuming the distance to be the 
maximum, 50 ft, the formula reads 

a = 473 u* 2 , (4.2-2) 

with a in pg and ut in deg/sec. Thus, a 0.2 deg/sec rotation 
rate (that planned for use during star acquisition for on-orbit 
gyro cal) would produce about 19 pg of error. This amount is 
not intolerable, but it can probably be avoided. 

Impulsive events , unless something is permanently 
ejected from the Orbiter, come in opposing pairs. Something 
begins to move, causing an impulsive event, and sometime later 
it stops, causing the opposing impulsive event. If both of 
these paired events occur during accelerometer calibration. 
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they will contribute no error. If accelerometer calibration 
begins or ends between the pair of events, then errors will 
occur. Again assuming a 200,000 lb Orbiter, an unpaired im- 
pulsive event during accelerometer cal will cause an error of 

6a = p/2059, (4.2-3) 

where 6a is in pg and p is the momentum of the impulsive event 
in lb- ft/sec. Obviously, it takes a substantial impulsive 
event to cause significant error in the accelerometer calibra- 
tion process; a 1 ton hatch moved at 5 ft/sec would produce 
only 5 pg of error. 

The only accelerometer instability worthy of note is 
a heading sensitivity error . The behavior of the accelerometer 
depends to some extent on the attitude of the cluster relative 
to the 1MU case. Ground calibrations maintain good repeatabil- 
ity by always calibrating with the same sequence of cluster-to- 
case attitudes. By and large, however, the heading sensitivity 
for the Orbitor 1MU is in the accelerometer scale factors, (see 
Chapter 2) and they are unexcited on orbit. Heading sensitivity 
of the biases is in the range of 7.5 pg rms, entirely tolerable 
for on-orb? t calibration. 

In summary, it appears that control of disturbing 
dynamical effects: 

• Thrust to an rms value of 2 lb or less 

• Vehicle rotation rate to less than 0.15 
deg/sec 

• Unpaired impulsive events to less than 
20,000 lb- ft/sec 

will allow accelerometer bias calibration on orbit to a level 
or roughly 20 pg, as Table 4.2-1 illustrates. 
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TABLE 4.2-1 

SUMMARY OF ON-ORBIT ACCELEROMETER 
BIAS CALIBRATION ERRORS 


ERROR SOURCE 

SOURCE MAGNITUDE 

RMS CONTRIBUTION TO 
CALIBRATION ERROR (pg) 

Quantization 

2 pulses 

2.55 

Thrust 

2 lbs rms 

10.00 

Rotation 

0.15 deg/sec 

10.64* 

Impulses 

20,000 lb- ft/sec 

9.71* 

Heading Sensitivity 

7.5 pg rms 

7.50 

RSS TOTAL 


19.24 


Peak values for worst geometry 


4 . 3 GYRO CALIBRATION 

In gyro calibration on orbit, a measured change in 
gyro drift angle divided by the time interval over which it 
accumulated is used to estimate the gyro drift rate bias 
(Ref. 3). To express the calibration in mathematical terms, 
assume the calibration period begins with time zero (t=0) and 
ends a period T later (t=T). In symbols, let: 

40 = measured change in gyro 

drift angle between t=0 
and t=T 

0(T) = estimated gyro bias drift 
rate at test’s end, t=T 

Then the estimation process is simply 
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0(T) = (4.3-1) 

To proceed, it is necessary to assume some statistical 
model of gyro drift rate. Ignoring for the moment heading 
sensitive drift rate effects, a model suitable for on orbit 
conditions (assuming no sensed acceleration) is: 

® U> = ^initial + rt + w(t> ’ (4.3-2) 

where b is the fixed bias drift rate at the initial time, r is 
a fixed ramp drift rate, and w(t) is a random walk process. It 
should be noted that the instantaneous drift rate bias for this 
model is 0(t) itself and the bias to be estimated is the final 
bias, 6(T). 


kfinal * ®< T > (4 ' 3 - 3 > 

The rms error in estimating the bias at the end of the test 
i s then 

oh 

o = E { (0(T) - 0(T) ] j (4.3-4) 

where E is the usual ensemble expectation operator. Without 
loss of generality, it may be assumed that th«. gyro drift angle 
is zero at the beginning of calibration (©(0) =0), so that 
the change in drift angle during the calibration period is just 
the final value 0(T). Then write: 

A0 m = 0 (T) + 6© m , (4.3-5) 

were 6© m is the measurement error in measuring the difference 
between drift angles at the beginning and end of the calibra- 
tion . 
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The differential equation describing the random walk 
process in Eq. A. 3-2 is: 

w(t) = n(t) , (A. 3-6) 

where n(t) is a Gaussian white noise process with covariance 
function 


E{n(t)n(i )} = o 2 6(t-t) (A. 3-7) 

where 6 is the Dirac delta "function." Also, without loss of 
generality, it will be assumed that the random walk process is 
zero at the beginning of the calibration interval (w(0) = 0). 

Making the reasonable assumption of statistical inde- 
pendence of the random walk white noise and of the gyro drift 
angle measurement error ( from each other and from the determin- 
istic parameters) yields for the mean square estimation error 

o 2 = <o m 2 + o n 2 T 3 /3 + r 2 T 4 /A)T 2 , (A. 3-8) 

2 

where a is the mean square angular error in measuring the 

ui 2 

drift angle change over the calibration interval, i.e., is 

the variance of 60 . 

m 

To evaluate the rms estimation error, a, and to opti- 
mize the calibration interval, T, requires numerical values 
for o^, a R , and r. From data provided by Rockwell Internation- 
al (Ref. 2A), the per-axis error in measuring 1MU alignment by 
observing two stars with the fixed star trackers is 160 sec 

rms. Since the difference of two such measurements is used in 

2 

the calibration to measure the change in drift angle, o m de- 
pends on the extent to which star tracker errors at two differ- 
ent times are correlated. Becarse the measurements may be 
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separated by several hours, a conservative assumption is that 
the errors are uncorrelated, and thus 

o m 2 = 2 (160 sec) 2 = 3.95xl0" 3 deg 2 (4.3-9) 

Singer Kearfott Division selects gyroscopes for Shuttle 
Orbiter use from a much larger population of gyroscopes by means 
of a 2 hr and a 17 hr test. The 17 hr test effectively insures 
that the total uncertainty due to the fixed ramp and random drifts 
satisfies* 


(119/13) a n 2 + Jl4. 75r| 2 < (0.013 deg/hr) 2 (4.3-10) 

Thus, firm limits on the two terms (each assuming that the 
other is zero) are: 


< (8.81xl0‘ 4 

deg/hr 2 ) 2 . 

(4.3-11) 

: < 1.85xl0' 5 

deg 2 /hr 3 . 

(4.3-12) 


The population of selected gyroscopes is likely to have flat 
distributions of ramp and white noise intensity paramters, so 
the mean square values of these parameters will be one third 
of the maximum square values. To split the error into the two 
pieces, consider the two extreme cases shown in Table 4.3-1. 
The rms gyro bias estimation error a (in deg/hr) may now be 
written in terms of the test interval T (in hr) directly from 
Eq. 4.3-8 as follows: 


Ramp Case 


o 


( 1/T) (3 . 95*10” 3 + 


6.47xl0‘ 8 



(4.3-13) 


Assuming the gyro drift rate model of Eq. 4.3-2 and the 17 hr 
test description of Ref. 25. 
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TABLE A. 3-1 

EXTREME CASES FOR DIVISION OF GYRO ERRORS 


MEAN SQUARE 



VALUE 

RAMP CASE 

RANDOM WALK CASE 

r 2 

2.59xl(T 7 deg 2 /hr 4 

0 

% 2 

0 

6.17xl0“ 6 deg 2 /hr 3 


Random Walk Case 


% 

a = (1/T) (3.95xl0“ 3 + 2.06xl(f 6 T 3 ) 


(4.3-14) 


The calibration interval may be optimized for either of these 
two cases. The results are: 

Ramp Case 

T opt = < 4 % 2 /r 2 > 1/4 = 15-7 hr (4.3-15) 

Random Walk Case 

T opt * <6a m 2 /°n 2 > 1/3 = 15 - 7 hr (4.3.16) 

The optimum test periods are identical in the two cases, to 
the significance shown. They are nearly equal because they 
happen to be close to the gyro test period of 17 hr used to 
allocate the errors to one case or the other. If the measure- 
ment uncertainty, a m , were greater, both optimum test periods 
would be longer than 15.7 hrs, but the test period for the 
random walk case would be the longer of the two. The opposite 
would happen for a smaller measurement uncertainty. 
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The sizes of the calibration errors in the two cases 
are not so nearly equal. At the optimum calibration time the 
rms gyro bias calibration errors are: 


Ramp Case 


a = 0.00565 deg/hr 


(4.3-17) 


Random Walk Case 

a = 0.00695 deg/hr (4.3-18) 

Thus, with a calibration interval of about 16 hrs, a one-sigma 

gyro bias calibration error should be less than 0.007 deg/hr. 

It should oe noted that the calibration accuracy is fairly 

insensitive to the calibration period. Assuming the error 

2 2 

model statistics (o , r, <i ) remain the same, the one-sigma 

m n 

error remains less than 0.01 deg/hr (in either case) for periods 
from 7 to 36 hrs. 


The foregoing analysis ignores the heading sensitive 
gyro drift effects that are potentially of large impact . Meas- 
urements on the prototype Shuttle 1MU showed rms variations in 
drift rate of 0.0120 deg/hr as the cluster-to-case attitude 
was varied 360 deg about the vertical gyro axis (Ref. 16). 

( ."her tests have produced drift rate changes as large as 0.07 
deg/hr with cluster-to-case attitude changes. Additionally, 
drift rate sensitivity to the direction of gimbal rotation and 
to the direction the gyro was last torqued for slewing have 
been noted. Software and hardware modifications to reduce the 
magnitude of these effects have been implemented, but a sizable 
residual heading sensitive drift will remain. 

Since the recommended length of on-orbit gyro calibra- 
tion (ignoring heading sensitivity errors) is about 10 orbital 
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periods, it might be expected that substantial averaging of 
drift rate variations due to cluster-to-case attitude sensitiv- 
ity would occur. The actual averaging obtained would depend 
on the details of Orbiter attitude during the calibration per- 
iod. If the Orbiter were inertially stabilized, cluster-to- 
case attitude would be constant, and no averaging would occur. 
Given the lack of hard information on expected magnitudes of 
heading sensitivity, and given the complex dependence of gyro 
drift rate on unpredictable cluster- to-case attitude, it is 
very difficult to assess the impact of heading sensitivity on 
gyro bias calibration accuracy on orbit . A rough estimate is 
that the accuracy could degrade to the 0.02 deg/hr (lo) range. 
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5 . ENTRY AND PRELAND NAVIGATION PERFORMANCE 


In Ref. 9, TASC presented results of a study of a 
candidate baseline navigation filter for the entry, preland, 
and landing phases of the first Orbital Flight Test Mission 
(OFT-1). The filter analyzed in that report is very pessimistic 
in its estimates of initial condition uncertainties and exter- 
nal sensor performance. A revised filter, provided by NASA 
(Ref. 26) with reworked estimates of initial errors and with 
less conservative assumptions on sensor performance, was speci- 
fied for comparison with the filter in Ref. 9. Performance of 
the revised filter is presented in this chapter. 


5.1 REVISED NAVIGATION SCENARIO 

The trajectory and measurement schedule for the OFT-1 
mission as described in Ref. 9, (see also Volume I, Section 
2.2.2) with only a change in the time of initiation of baro- 
altimeter measurement processing, is assumed for this study 
(see Table 5.1-1). Previously, the baro-altimeter was switched 
on at an altitude of 84,000 ft (corresponding to a vehicle 
speed of Mach 2.2) and temporarily switched off when the vehicle 
experienced speeds between Mach 1.3 and Mach 0.9. This restric- 
tion prevents the processing of this altitude measurement until 
after the 31,000 ft altitude, and yields a total baro-altimeter 
operating time of 42 sec as compared to 183 sec for the previous 
filter. 
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TABLE 5.1-1 

KEY EVENTS FOR OFT-1 MISSION 




NAVIGATION 



ELAPSED 

TIME TO 

ALTITUDE 

RELATIVE 

RANGB-TO- 

PHASE 


EVENT 

TIME 

TOUCHDOWN 

ABOVE GEOID 

VELOCITY 

GO 



(MC) 

(*ec) 

(ft) 

(ip«) 

(ft) 

ENTRY 

r 

Entry Interface. T n 

0.0 

1630.1 

398,837. 

24296. 

19.600.000. 

NAVIGATION ’ 
PRASE 


u 

First Drag Measurement. T DRAQ 

168.2 

1441.9 

988,174. 

24399. 

15.400.000. 



First TACAN Measurement. T TACAN 

1062.2 

677.9 

144,780. 

7391. 

1.426.000. 



Last Drag Measurement 

184ft. 0 

382.1 

100.263. 

3446. 

478,000. 

PRELAND 
NAVIGATION « 


First Baro Altimeter Measurement. 
t ba 

First MLS Azimuth, Elevation, ft IMS 

1493.6 

136.3 

31.230. 

745. 

67.268. 

PHASE 


1636.0 

94.1 

30.370. 

662. 

51.000. 

Measurements. Tuie 
(No TACAN Measurement) 


w 

(No Baro Altimeter) 






, 


Reduce Filter Dimension. Tgy 

1662.9 

67.2 

13.680. 

622. 

35.000. 

LANDING 
NAVIGATION ^ 


Switch from 3.64 sec Update Cycle 
to 1.92 sec Update Cycle 







First Radar Altimeter Measurement. 

1622.4 

7.7 

2.387. 

382. 

2.700. 

ilAf AHMI AWil 

PHASE 


Tra (No MLS Elevation Measurement) 
Switch from 1.92 sec Update Cycle 








to 0.16 sec Update Cycle 








Touchdown 

1630.1 

0.0 

2.293. 

315. 

0. 


The initial filter covariance matrix has been reworked 
by NASA and may more closely reflect the actual navigation 
errors at deorbit burn. The previous filter assumed initial 
position and velocity errors of equal magnitude in each com- 
ponent of position and velocity. However, some components are 
more difficult to track during orbital operations than others 
(e.g., downrange position is more difficult to estimate than 
crossrange position) and the revised filter covariance takes 
such differences into account. In addition, the initial mis- 
alignments specified for the revised filter are somewhat smaller 
(approximately 0.2 deg) th*n those used previously (0.25 deg). 
Table 5.1-2 presents the initial condition uncertainties in 
the radial, downrange, and crossrange (R,DR,CR) directions. 

The initial covariance matrix also assumes a -0.9 correlation 
between downrange position and vertical velocity. 
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TABLE 5.1-2 

INITIAL CONDITION ERRORS FOR REVISED FILTER 

T-0811 


STATE 

R 

DR 

CR 

Position (ft) 

2500 

25000 

2500 

Velocity (fps) 

26.3 

3.3 

3.3 

Misalignments (mrad) 

3.464 

3.464 

3.464 


The navigation filter for the entry through landing 
phases is a variable dimension Kalman filter with up to 12 
states. After Tgy (approximately 12,000 ft altitude) it is a 
6 state complementary filter. A schedule of filter states is 
presented in Fig. 5.1-1. 


APPROXIMATE 

,45< j» Fr 


100.000 FT 


31,000 FT 70.000 FT 

I I 


12,000 FT 100 FT 

I I 


R 11206b 
TOUCHDOWN 

I 


filter 
STATES 


16 


k 


12 


k 


POSITION AND VELOCITY 


MISALIGNMENTS 


. DRAG UPDATE , 

10 h — H 


BARO 

altimeter 

i I MLS ELEVATION . 

N — »4«— — M 


TACAN DM6 


^ H 


TACAN VOB 


, MLS AZIMUTH , 

-k« — — H 


Figure 5.1-1 Schedule for Entry, Preland, and Landing 

Phase Navigation Filter States 


r .-3 
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The gain computation algorithm in the Kalman filter 
mode (i.e., prior to Tgy) includes measurement underweighting 
in order to minimize the possibility of filter divergence: 


K 


PH 


u 7Q HPH T + R 


-1 


(5.1-1) 


where the measurement underweighting factor u^q is 1.2 when 
the rss (filter-indicated) position error is greater than 3280 
ft and 1.0 otherwise. 


The sensor parameter values for the entry and preland 
navigation phase Kalman filter are given in Table 5.1-3. The 
major changes from the filter in Ref. 9 include: 

1) A reduction in the rms drag-altitude pseudo- 
measurement noise from 5148 ft to 2200 ft, 
and 

2) Greatly improved TACAN bearing statistics. 

The standard deviation of the TACAN bearing correlated error 
has been decreased from 29 mrad to 9.64 mrad, and the measure- 
ment noise has been halved to 7 mrad. 


The baro-altimeter error statistics for the filter 
are poorer than those used previously, but the baro-altimeter 
is operated for a much shorter time, thereby diminshing the 
effects of the poorer performance estimates. The TACAN range 
parameters are slightly better than, and all the MLS error 
statistics are exactly the same as, those used in the previous 
filter. 


The time history of the filter-indicated performance 
between deorbit burn and initiation of the landing phase is 
summarized in Table 5.1-4. These results represent what 
optimal performance would be if this 1 state filter model 
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TABLE 5.1-3 

FILTER STATE AND MEASUREMENT ERROR STATISTICS FOR 
ENTRY AND PRELAND NAVIGATION FILTER 


T-0812 


MEASUREMENT 

CORRELATED ERROR FILTER STATE 

MEASUREMENT 

NOISE 

STANDARD 

DEVIATION 

CORRELATION 

TIME 

(SEC) 

STANDARD 

DEVIATION 

Drag-Altitude 

9842 ft 

800 

2200 ft 

Baro Altimeter 

2462 ft 

400 

868 ft 

TACAN DUE 

444 ft 

400 

133 ft 

TACAN VOR 

9.64 mrad 

400 

7 mrad 

■ 

80 ft 

0» 

26 ft 


0.96 mrad 

• 

0.3 mrad 


0.86 mrad 

«r 

0.3 mrad 


TABLE 5.1-4 

FILTER- INDICATED PERFORMANCE 


T- 04 94-1 


EVENT 

ELAPSED 

TIME* 

RMS POSITION ERROR 
(ft) 

RMS 

VELOCITY 

(fps) 

ERROR 


(sec) 

R 

DR 

CR 

R 

DR 

CR 

Initiation of Deorbit 
Burn 

— 

2500 

25000 

2500 

26.3 


3.3 

3.3 

Entry Interface 

0.0 

16505 

24859 

8360 

38.5 


18.5 

12.0 


163.2 

19540 

27289 

9628 

m 

1 

20.2 


First Drag Measurement 

188. 2 * 

11577 

26128 

.9809 

m 


16.5 



1029. 1* 

9792 

29732 

28194 

m 

1 

37.9 

M 

First Tacan Measurements 

1052. 2 + 

9484 

12977 

17373 

EB 

1 

32.7 

SB 


1489. 9 + 

1966 

1182 

991 

14.2 


8.7 

11.2 

First Baro Altimeter 
Measurement 

1493.8* 

1472 

901 

822 

13.9 


8.1 

10.4 


1536.0" 

732 

Bi 

262 

11.1 


8.2 

8.0 

MLS Acquisition 

1536.0* 

57 

■9 

65 

11.0 

J 

6.8 

5.6 

Initiation of Landing 
Navigation Phase 

150? 9" 

44 

78 

45 

2.2 

2.5 

2.3 


’Add 1187.0 sec for elapsed time from deorbit. 
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were a complete and accurate representation of real-world error 
dynamics. The major difference from the behavior observed in 
Ref. 9 is a marked improvement in all components of the pre- 
dicted position and velocity errors through the period of TACAN- 
only operation (i.e., 1052.2 sec to 1493.8 sec). This improve- 
ment is expected since the TACAN error statistics, particularly 
bearing, are improved over the previous filter values. At the 
switch to MLS, however, velocity errors are larger than those 
previously indicated because of the baro-altimeter ' s shorter 
operating time and poorer measurement statistics. At the switch 
to the landing navigation phase these velocity errors are 
approximately twice those indicated in Ref. 9. 


5.2 NAVIGATION SYSTEM PERFORMANCE 


Performance results for the Space Shuttle entry and 
preland navigation system are summarized in this section in 
the form of a table showing rms position and velocity errors 
at specific flight times, and in the form of partial error 
budgets at the end of the entry phase ( ) and at the 
initiation of MLS measurement processing < s > • These 
tables are based on only those error sources which have been 
major contributors to navigation performance down to the 
29,000 ft altitude in recent TASC studies . In particular these 
error sources are 


• Initial position and velocity errors 
(Group 1) 

• Initial platform misalignments (Group 1) 

• Drag-altitude, and TACAN range and bearing 
measurement noises (Group 1) 

• Accelerometer biases (Group 2) 

• TACAN range bias (Group 10) and bearing 
bias (Group 19) 


5-6 




THE ANALYTIC SCIENCES CORPORATION 


• Baro-altimeter bias, scale factor and 
static defect (Group 12) 

• The 1962 standard atmosphere modeling 
error (Group 21) 

The error magnitudes associated with these error sources are 
the same as those used in Ref. 9. While several of the addi- 
tional error sources could become important for a particulai 
filter mechanization , it is reasonable to assume that they 
would add only a few percent to the overall performance esti- 
mates. 


Table 5.2-1 presents a time history of position and 
velocity errors for the revised navigation filter (generated 
by root -sum- squaring the individual contributions of the error 
sources listed above). The principal difference from the per- 
formance shown in Ref. 9 is an increase in vertical and down - 
range errors below an altitude of 85,000 ft, presumably because 
baro-altimeter measurements are suppressed until late in the 
mission . The initial covariance for the revised entry filter 


TABLE 5.2-1 

ESTIMATED PERFORMANCE FOR REVISED NAVIGATION FILTER 

1 = 8812 . 


EVENT 

IT TTTTTnr 

ELAPSED 

•mir* 

RMS POSITION ERROR (ft) | 

«” tEugm — »»> i 

Avl A i V& /£, 

(ft) 

(sec) 

R 

DR 

CR 

R 

DR 

CR 

Before First 
Drug ileas. 

287,000 

188.2* 

2110. 

11520. 

1366. 

14.36 

2.26 

1.69 

Before First 
Tacar. Xe&s. 

145,000 

1052.2* 

2626. 

12391. 

5076. 

19.64 

9.45 

13.97 

After First 
Tacan ilea s. 

145,000 

1052.2* 

2790- 

4279. 

14918. 

16.61 

13.44 

29.72 

After Last 
Brag 

1C?, 000 

1248.0* 

m 

874. 

3159 

15.86 

10.54 

22.83 

After First 
Bars We a s. 

31,000 

1493. 8 r 

1865. 

835. 

1078. 

11.81 

7.38 

8.84 

Before First 
ULS lie as. 

20,000 

1536.0* 

829. 

807. 

366. 

11.93 

9.64 

5.99 


"Add 1X87.0 sec for Elapsed Tine fros Deorbit. 
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is still pessimistic relative to the truth model covariance as 
is indicated by a comparison of Table 5.2-1 with the filter- 
indicated performance. Table 5.1-4. 

The initial pessimism of the revised filter along 
with the lower values of predicted drag-altitude pseudo-measure- 
ment noise yields larger radial position errors during blackout 
than were experienced in Ref. 9. This behavior suggests that, 
for the assumed truth model error magnitudes, the revised filter 
is overweigting these altitude measurements. Similarly, the 
immediate increase in crossrange errors with the processing of 
the first TACAN bearing measurement is due to the fact that 
the TACAN bearing error statistics in the filter are now smaller 
than those in the truth model , again leading to possible over- 
weighting of the bearing measurements. Finally, as discussed 
in Ref. 9. the large velocity errors are primarily due to the 
large 1MU misalignments. 

Table 5.2-2 presents a partial error budget for the 
entry navigation phase. The only major contributor (as deter- 
mined in Ref. 9) which is not included is the nonstandard aero- 
dynamics time-varying bias. This error source was a major 
contributor to the vertical position error, but its importance 
for this filter would be overshadowed by the contribution from 
the drag-altitude pseudo-measurement noise. Of the total pro- 
jected performance values only the vertical position error is 
considerably different (in this case, worse) than that of the 
previous filter. 

The 5000 ft standard deviation for drag-altitude 
pseudo-measurement noise used in Ref. 9 corresponds to a 20% 
random atmospheric density error. Although this number is 
consistent with values assumed in recent studies (Refs. 27 and 
28). inspection of available atmospheric density data (Ref. 29) 
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TABLE 5.2-2 

REVISED PARTIAL ERROR BUDGET FOR 
THE ENTRY NAVIGATION PHASE 


T-0814 





RMS NAVIGATION ERROR AT 145.000 ft | 


ERROR SOURCE 

VALUE 

POSITION (ft) 1 

VELOCITY (fps) | 




R 

DR _ 

^CR 

^ R 

DR 

CR 

2. 

FILTER STATES AND MEASUREMENT 
NOISE 

Initial Condition Errors 
Positions and Velocities 

R*f. 9 

224. 

11100. 

3752. 

14.75 

3.25 

2.17 


Misalignments 

265 sec 

ISO. 

970. 

2416. 

6.07 

6.46 

13.54 


Uocorrelated Measurement Noise 
Altitude 

5000 ft 

2603. 

2037. 

1257. 

10.61 

4. 76 

1.64 


INS .Quantization Noise 

(1250 ft) 
0.0326 fps 

(651. 

^ -r 

} (509.) 

(314.) 

(2.65)(1 .19)(0.41 1 

11. 

I MU- RELATED STATES 

Accelerometer Errors 
2. biases 

50 tig 

220. 

1903. 

1617. 

1.52 

1.66 

1.22 


3. Scale Factors 

100 ppm 



- 

- 

„ 



4. asymmetries 

100 ppm 
IS deb 

— 

— * 

- 





5. Nonorthogonalities 



_ 



_ 


Gyro Errors 

7. Biases Drifts 

0.035° /hr 




i 




6. Mass Unbalances 

0.025°/br/g 

— 




— 

— 

IV. 

DRAG- RELATED STATES 
21. Non-Standard Deasity 

1962 Standard Atmosphere 

4-TeraC <9%) 

5. 

4649. 

1266. 

4.06 

3.40 

1.75 


Modeling Error 
Time- Varying Bias 

Ref. 9 
Vinter (<8l) 

— 

— 

— 

— 

— 

— 


First-Order Markov 

( <6%) Ref. 9 




_ 


— 


22. Mon-Standard Wind 
Westerly 

Vinter 








Crosswind) 

(<295 fps)Ref .9 
(<110 fps) 








Headwind f 
23. Non-Standard 
Aerodynamics 

Time- Vary in; Bias 

Ref. 9 

i 

«9.5*) 



i 

— 





First-Order Markov 

Ref. 9 
51 of C^ 

— 

— 

i 

— 

— 

— 

Total Projected Performance 


2626. 

12391. 

5076. 

19.64 

9.45 

13.97 

Tcval Projected Performance With 
Alternate Measurement Noise 


(?38.>a2233. >(4928. ) 

ITir7TT7w^Tt^y^y 


•Negligible contributor 
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does not support a value greater than 1250 ft (5%). Because 
the revised filter is much more sensitive to the magnitude of 
the drag-altitude pseudo-measurement noise than previous fil- 
ters, it is desirable to include in Table 4.2-2 the more opti- 
mistic (and, perhaps, realistic) vertical channel performance 
which would be attained with this smaller value. The parenthet- 
ical entries in the error budget provide both error contribu- 
tions and total projected performance with the smaller measure- 
ment noise. Even with the lower value of measurement noise, 
this error source is at least a major contributor to the verti- 
cal position error. There is a clear need for adequate test 
data and a better understanding of this error source . 

Except for the contribution of the initial misalign- 
ments to the vertical position error, which is much smaller 
than previously observed, all the remaining contributions listed 
in Table 5.2-2 are approximately equal to those of the previous 
filter. 


A partial error budget for the revised navigation 
filler at 20,000 ft is given in Table 5.2-3. Because of a 
lack of reliable data on TACAN bearing errors, alternate values 
of the TACAN bearing errors, and the corresponding total pro- 
jected performance using these alternate values, are also 
included in the table. The total projected performance in 
nearly all channels is somewhat poorer than that of the pre- 
vious filter. 

The filter's heavy weighting of the TACAN bearing 
measurement is indicated by the very large contributions of 
the TACAN bearing bias and measurement noise to all components 
of both the position and velocity errors. Even with these 
error sources set at the lesser values, i.e., 6 mrad bias and 
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TABLE 5.2-3 

REVISED PARTIAL ERROR BUDGET AT 20,000 FT 


T-0S03- 1 




RMS NAVIGATION ERROR AT 20,000 


ERROR SOURCE 

VALUE 

POSITION 

(tt> 

VELOCITY ( fps) ] 




t m DR 

CR 


on 

_ cn 

I. Fll/TEH STATES AND MEASUREMENT 

KSTSE 

Initial Conditions 








Positions and Velocities 

Ref. 9 

29.1 

31.5 

13.9 

4.04 

0.25 

1.02 

Misalignments 

265 sec 

77.4 

103.2 

49.7 

3.89 

1.73 

3.93 

Uncorrelated Measurement Noise 








TACAN Range 

100 ft 

141.2 

62.1 

31.7 

1.93 

0.57 

0.54 

TACAN Bearing 

12 mrad 

348.7 

429.8 

133.4 

4.61 

5.50 

1.49 


(6 mrad) 

0.24 e h ' 24018 ft 

tl74.4)<214.9) 

(66.7) 

|(2.30)(2.75) 

(0.74) 

Baro- Altimeter 

if 

_ 

_ 




_ , , 

INS Quantization Noise 

0.0326 fps 

— 

— 

— 

— 

— 

— 

11. I MU-RELATED STATES 
Accelerometer 








2. Biases 

SO lag 

9.5 

12.6 

6.0 

0.65 

0.19 

0.52 

3e Scale Factors 

100 ppm 

— . 


— 

_ 


__ 

4. Asymmetries 

100 ppm 
IS £eb 





_ . 

- 

5. Honor tbogonali ties 




— r 

— n» 


Gyro 








7. Bias Drifts 

0.035° /hr 




— - 



8. Mass Unbalances 

0.02S°/br/g 

— 

— 



— 

— 

IV. EXTERNAL AID-RELATED STATES 
10. TACAN Range Bias 

385 ft 

449.4 

150.6 

70.4 

4.96 

1.90 

2.61 

11. TACAN Range Scale Factor 

1*00 «ppm 







12. Baro-Al time ter Errors 








Bias 

Scale Factor 
Static Defect 

100 ft ) 

3% 4 ( 

1.52 * 10'* ( 
f t/f t 2 /sec 2 ) 
20 ft 

69,7 

34.5 

20.4 

J.D5 

0.30 

0.50 

First-Order Markov 







16. TACAN 8urvey Errors 

1 ft 



. 

— 



19. TACAN Bearing Bias 

12 mrad 

575.5 

652.6 

326.9 

7.74 

7.46 

3.10 


(6 mrad) 

(267. 7) (326. 4) 

(163.4) 

3.87X3.73) 

(1.55) 

1. ,21., 22. ,23. 








Drag-Related Errors 

Ref. ? 

— 

— 

— 

— 

— 

— ~ ' 

Total Projected Performance 


628.6 

806.5 

365.7 

LI. 93 

9.64 

5.99 

Total Projected Performance With 
Alternate TACAN Bearing Model 


(589.0X438.5) (200.5) 

;e.02)(5.35)(5.20) 


'Negligible contributor 
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6 mrad measurement noise, they are major contributors to all 
position and velocity errors. As is the case with the alti- 
tude pseudo-measurement noise, the importance of the TACAN 
bearing errors in this filter underlines the need for a better 
understanding of the true TACAN bearing error mechanisms and 
their magnitudes . 

The TACAN range errors are major contributors to verti- 
cal channel errors. These contributions are somewhat larger 
than observed for the previous filter, due to the slightly im- 
proved TACAN range error statistics and the delayed utilization 
of the baro-altimeter . The baro-altimeter error sources have 
been reduced to insignificant contributors because the baro- 
altimeter is operated for a shorter length of time and the 
filter has much poorer error statistics than the previous fil- 
ter. 


A principal difference between these results and those 
of the previous filter is the contributions of the initial 
position and velocity uncertainties. Where previously this 
error source was a major contributor to all position and veloc- 
ity errors, it is now a major contributor to only the vertical 
velocity error. The revised filter, with the reworked estimates 
of the initial conditions, does not generate the large misalign- 
ment estimation errors seen in Ref. 9. As a result, the mechan- 
ism which couples these initial position and velocity uncertain- 
ties to the position and velocity errors, through the misalign- 
ments, does not produce the substantial error contributions 
observed in Ref. 9. 

The initial misalignments continue to be major con- 
tributors to vertical and crossrange velocity errors. This is 
expected since the initial misalignment estimates were reduced, 
but are still quite large. 
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Although error budget runs for all the error sources 
considered above were made down to touchdown, the corresponding 
error budget tables are not included here. The changes in 
these tables from the results observed in Ref. 9 are minimal. 

The principal change is that the TACAN bearing errors are major 
contributors to crossrange position and velocity errors through- 
out the remainder of the flight. 

Overall the performance of the revised filter down to 
the 20,000 ft altitude is somewhat poorer than that indicated 
for the previous filter, primarily because baro-altimeter meas- 
urements are suppressed until late in the mission. At altitudes 
below 20,000 ft, where MLS- and radar altimeter-related errors 
dominate, the performance of the two filters would be approxi- 
mately the same. The revised filter experiences two important 
differences from the previous baseline filter: 

• The altitude pseudo-measurement and the 
TACAN bearing measurement assume very 
significant roles during the entry and 
preland navigation phases. 

• The initial position and velocity un- 
certainties do not yield the large mis- 
alignment estimates which plagued the 
filter in Ref. 9. 

It should be further noted that the altitude pseudo-measurement 
noise and the TACAN bearing errors are not well documented or 
well understood. Their importance in this candidate filter 
und erlines the need for improved (better unde rstood) error 
models for these error sources. 
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6 . 


SUMMARY AND CONCLUSIONS 


This volume of this final report presents the results 
of four distinct study areas in Space Shuttle navigation asso- 
ciated with the baseline navigation system configuration. A 
brief summary of the major conclusions arrived at in the pro- 
ceeding four chapters is presented here. 


6.1 1MU CALIBRATION AND ALIGNMENT COVARIANCE ANALYSIS 


Covariance analysis techniques were applied to develop 
a detailed error budget for the 0FT/1MU calibration and align- 
ment algorithm. Similar results were presented in Refs. 4 and 
5, but the present results apply to a revised version of the 
algorithm (Ref. 3, Dec. 1976) with 13-position Hangar Cal A 
and 7-position Preflight Cal A sequences designed to avoid 
transient effects of accelerometer gain switching. Several 
other changes in procedures and analysis, as are detailed in 
Chapter 2, are also incorporated. 

Complete error budgets for the Hangar and Preflight 
calibrations are presented. Significant conclusions to be 
drawn from the error budgets include: 


• The specified Hangar Cal A calibration 

accuracy appears to be achievable, except 
for low gain accelerometer bias and accel- 
erometer scale factor asymmetry. Acceler- 
ometer scale factor errors sensitive to 
gimbal attitude are primarily responsible 
for the lack of total compliance. 
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• The specified Hangar Cal B calibration 
accuracy appears to be achievable. 

• The specified Preflight Coarse Alignment 
accuracy appears achievable. 

• The specified Pre flight Cal A accuracy 
appears achievable, except for low and 
high gain accelerometer bias. Again, 
the primary reason for noncompliance 
lies with the accelerometer scale factor 
attitude sensitivity. 

• The specified Gyrocompass Goodness Test 
accuracy appears achievable only if the 
gyro attitude sensitive drift can be 
essentially perfectly compensated. 

• The specified Gyrocompass Azimuth Align- 
ment accuracy appears achievable only if 
the gyro attitude sensitive drift is 
almost perfectly compensated. 

• The specified Velocity and Tilt Initiali- 
zation accuracy appears achievable. 


A more detailed presentation of conclusions and recom- 
mendations for 1MU calibration and alignment are to be found 
within Chapter 2 of this Volume. The principal recommendations 
are for further validation and improvement of the IMU sensor 
error model terms that heavily impact both the calibration and 
alignment procedures and the subsequent Orbiter navigation 
performance. 


6.2 POSTFLIGHT IMU ERROR RECOVERY FOR APPROACH AND 
LANDING MISSION PHASES 


Atmospheric flight tests such as ALT provide an oppor- 
tunity to evaluate both the navigation performance achievable 
with the IMU and the true calibration and alignment uncertain- 
ties in a flight environment. Postf light recovery of inflight 
IMU errors, especially platform alignment errors, would provide 


6-2 


THE ANALYTIC SCIENCES CORPORATION 


an accurate reference for the evaluation of other Shuttle sub- 
systems, such as the air data system. The study of postflight 
error recovery focuses on the test design, on the test instru- 
mentation, and on the data analysis techniques involved. A 
computer simulation was created to examine these issues. 

The results of the simulation indicate a potential 
for a significant reduction in the preflight uncertainties of 
level axis errors (tilts, level drifts, and accelerometer 
biases). Azimuth axis uncertainties are also reduced, but not 
so significantly. Certain errors, including accelerometer 
scale factors, are essentially unobservable, while others, in- 
cluding gyro bias and input axis mass unbalance for the azimuth 
axis, are inseparable. Nonetheless, it appears that postflight 
processing of ALT measurements could maintain IMU alignment 
errors at a level smaller than the preflight misalignment values . 


6.3 ON-ORBIT CALIBRATION OF SHUTTLE IMU ACCELEROMETER 
AND GYROSCOPE BIASES 


Accelerometer biases and gyroscope bias drifts may be 
calibrated in the sensed-acceleration free environment on orbit. 
A short (320 sec) accelerometer calibration will be accurate 
to about 20 Mg (lo) , provided that certain acceleration pro- 
ducing events (venting, thrusting, rotation, and impulsive 
events) can be held to what appear to be easily achievable 
levels. Accelerometer stability itself is of little concern 
over so short a calibration interval. 

On orbit gyroscope calibration, on the other hand, is 
largely a tradeoff of instrument stability and startracker 
measurement error, the startrackers being used as the attitude 
reference for gyroscope calibration. With some uncertainty 
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engendered by the potential effects of heading sensitive gyro 
drift, it appears that a 16 hr gyro calibration period would 
achieve bias drift calibrations to a level of 0.007 deg/hr 
(lo) . Acceleration effects important in on-orbit accelerom- 
eter calibration have no relevance to on-orbit gyro calibration. 


6.4 ENTRY AND PRELAND NAVIGATION PERFORMANCE 

This study is a revision of previous work (Ref. ?) 
with some changes in assumptions. The main issue is naviga- 
tion filter performance in the entry, preland, and landing 
phases of the first Orbital Flight Test mission. The major 
changes in assumptions from the earlier study are better drag 
updating and improved TACAN performance in the present study, 
but a lowering in the maximum altitude at which the baro alti- 
meter is assumed usable (associated with passage of the Orbiter 
through Mach 1 ) . 

The present results show an increase in vertical and 
downrange errors below a 85,000 ft. altitude because no baro 
altimeter measurements are yet available , and a greater de- 
pendence on the accuracy of the drag updates. This degraded 
behavior of the filter continues even down to 2000 ft (the 
baro altimeter being assumed available at 31000 ft), showing 
the importance of the early altitude measurements . Because 
the drag update errors and the TACAN errors are not well docu- 
mented, and because this revised filter depends so strongly on 
these errors, a better understanding of the statistical behav- 


ior of drag updating and TACAN bearing errors should be pursued . 
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